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Introduction. 


FeséR has shown ¢ that the Fourier development corresponding to a function 

J(«), of period 27, will be summable ¢ to the value 
+ 0) + f(x — 9)] 

at all points at which the function f(a) is continuous or has a finite jump, if 
J(«) is finite and integrable or even if f(x) becomes infinite at a finite number 
of points, provided it remains integrable. He has further shown in the article 
referred to that the development is uniformly summable throughout any closed 
interval which does not include a point of discontinuity of the function. He 
has used these properties in connection with a general theorem about conver- 
gence factors to establish the fact that the formal results obtained in the discus- 
sion of many problems in Mathematical Physics which involve the development 
of an arbitrary function in a Fourier’s series really furnish a solution of the 
problem, even when the development of the function is a divergent series. 

FEJ&éR’s idea of investigating the nature of the divergence of a development 
of an arbitrary function in terms of normal functions and then applying general 
theorems about convergence factors to determine the behavior of the series when 
such factors are introduced, can be applied to many other developments that 
occur in Mathematical Physics. The present paper is devoted to a study, from 
this point of view, of the developments in terms of Bessel functions. 

Although the formal work of obtaining the development of an arbitrary fune- 
tion in terms of Bessel functions goes back to FourIER, the rigorous discussion 


of the conditions under which the development is convergent has not made 


equal progress with that discussion in the case of the ordinary Fourier’s series. 
1 
This is doubtless due to the less elementary character of the functions involved 


* Presented to the Society, September 11, 1908. 
tCf. Mathematische Annalen, vol. 58 (1903-04), p. 51. 
¢ The word summable will be used in this article in a somewhat special sense. Cf. the defi- 
nition in ¢ 2. 
Trans. Am. Math. Soc. 26 391 
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and to the fact that the formule in the case of the Bessel functions are con- 
siderably more cumbersome. As far as I am able to learn the first and almost 
the only satisfactory investigation of the convergence of the development to the 
value of the function developed is due to Dint.* KNESER has discussed,t from 
the standpoint of integral equations, the convergence of the development in 
terms of J,, and has stated that similar methods will apply to the developments 
in terms of J, where vy> 90. A number of other discussions of the convergence 
of the development have been published, but I do not know of any treatment, 
besides those mentioned, in which the reasoning is complete. 

The discussion of the summability of the development, given in this paper, is 
analogous in its general outlines to the discussion of the convergence given 
by Harnack.{ Harnack’s treatment is not adequate but his method is 


suggestive. 


§1. Distribution of the roots of an equation involving Bessel functions. 


We wish to discuss in this section the distribution of the positive roots of the 


equation 


(1) + hd, (2) = 0, 


where /,() represents a Bessel function of the vth order, v is a constant posi- 
tive or zero, and / and / are any constants not both zero. 
We find it convenient to consider first the following trigonometric equation 


k 
(2) cos (A — a) +, sin(A—a)=0 ( k, a constants ) 


whose roots will be found to approximate to the roots of equation (1) for large 
values of 2. 

Lemma 1. Jf we represent by r,, the successive positive roots 
of equation (2), we have 


(m=1, 2, 3, ---) 


m 


(3) — 


where K, is a positive constant. 
This lemma is obviously true if £ = 0. Let us then first consider the case 


in which & > 0, and let A, be any positive value of A such that 
(4) cos (A, —a) = 0, sin(A, —a)=—1, 


(5) A, >&k. 


* Serie di Fourier (Pisa, 1880), pp. 246-269. 
t Mathematische Annalen, vol. 63 (1906-07), p. 505. 
t Mathematische Annalen, vol. 35 (1890), p. 41. 
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Since, in view of (4), cos(A— a) steadily increases from 0 to }V 3 and 
sin (A — a) steadily increases from — 1 to — } while A increases from A, to 
A, + 47, it follows that the expression on the left hand side of equation (2), 
which is negative for \ = A, and, in view of (5), is positive for A = A, + 47, 
must vanish at one and only one point in the interval A,=A=A,+ 47. 
Moreover, this expression is always positive in the interval A, + 47 =A=A,+7, 
and hence equation (2) has one and only one root in the interval A, =A=A,+7. 
If we denote this root by X’ , we have 


cos +7—a)+ 5, (A, — 


m 


(6) = sin(A, —a)+ sin(A, 
wk sin — a) 
rN (A +7) 


m m 


< 


Since cos (A — a) + (k/2) sin(A— is positive for + 77, it follows 
from (6) that ,, lies in the interval A, + 7<A<A, +7. Moreover we 
have from (6) 


wk sin + — a) 


+ 7) 
= cos(A, ,—a)+ sin(A; ,—@)—cos (A; +7—a)— 4+ (A, +7—a) 
(7) > cos — @) — cos (A, 4+ x7 Sin (A, + 7 — @) 
m+1 

_ 7 (Ai 

+r—X. 
= — cos( A) — ksin(X; 


But since from the Law of the Mean 


cos (A; ,, —a) —cos(A) + — | 
-a<&<4 
we have from (7) 
ksin(X\' +r7r—a r’ 
sin (A, + (A, + T ) 


and since, moreover, 


sin sin (A, + m7—a) = —sin(A) —a)> — sin(a, + a) =5 
and A; ., — A, — 7 is negative, we have finally 
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A 

4rk 

<2k| 

A A 


m m m 


m 


The inequality (8) has been derived under the assumption that ,, the root of 
cos (A — a) = 0 that immediately precedes 2X’, is such that sin (A, — a) = —1 
and the inequality (5) is satisfied. The first restriction was made for definite- 
ness and the proof of formula (8) can be carried through in a similar manner if 
(4) is replaced by 
(4’) cos (A, —a)=0, sin(A,—a)=1. 


Moreover, since there are only a finite number of positive roots of equation (2) 
lying in the interval 0 << X, where X, is the least value of A, satisfying (4) 
and (5), or (4) and (5), we may choose K, > 47k and such that the inequality 
(3) is satisfied for all positive roots of (2) that are less than r,- Hence, in 
view of (8), we have 


K 


and the lemma is proved for the case where k > 0. 

The proof for the case where / <0 is quite analogous, the principal differ- 
ence being that the successive roots of (2) come before instead of after the cor- 
responding roots of cos (A — a) = 0 to which they approximate. 

We will consider next the equation 


k 
(9) cos (A — a) +, sin(A—a) + €(A) = 9 


where & and a are constants and €() is a continuous function of 2 satisfying 


the inequality 
M 
(10) e(A)| < (M a positive constant). 
Lemma 2. Corresponding to each root X\, of equation (2) there is a root X, 
of equation (9) such that 


where K, is a positive constant. 
Let us first consider the case where 4 > 0; and let A, be any positive value 
of X greater than ?7 and such that 


(12) cos (A, —a) = 0, sin(A,—a)= —1; 
k 1 M T 


0 
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Then equation (2) will have one and only one root in the interval 
Ay <A<A,+ 447. Denote that root by»). For any positive increment AX 


we have 
h 
cos (A) + AX — a) + x + + Ar—a) 
(14) 
s(n’ + AX 
=cos(A, + — a) — cos ( — a) 
+> [sin (A’ + AX — a) —sin(A’) —a)}]. 
But 


cos + AX — a) — cos (A) — a) =— + PAX—a) (0<4<1), 


or if we take AX < }z so that 


+ — a) < sin 3 
Ar 
(14a) cos (A) + AX — a) — cos(A) —a)>-5 
Moreover, in view of (13), 
k Ar 
(14d) + AR) sin(A, +AX—a) < Ax < 
and since 
sin + AN— a) —sin (A, — a) = Adeos(A, + 6,AX—a) (0<4,<1), 


we have from (13) 
(14c) x7 {sin (A, + AA — — sin — @)| < 


m 


Combining (14) with (14a), (145), and (14c), we obtain 


+ AN—a) gy sin(X, + AN—a)>"y 


If now we choose AX so that 


(16) AA = 2< 49° 


we have, since + AX> A, — 
Ar M M 
= 
40 (A, + Ad)? 
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and hence, in view of (10) and (15), 


In a precisely similar manner it may be shown that 


k 
(18) — AX — a) + x’ sin (A, — Ar —a)+e(A 


A 


It follows from (17) and (18) that equation (9) has at least one root in the 
interval 
ARK + AX. 


If we represent such a root by \,, we have from (16) 


(19) 4M 1 9M 


The inequality (19) has been derived under the assumption that A,, the root 
of cos (A — a) = 0 that immediately precedes X’ , is such that sin (A, — a) = — 1 
and the inequalities (13) are satisfied. The first restriction was made for 
definiteness and the proof can be carried through in a similar manner if (12) is 
replaced by 
(12’) cos (A, —a) = 90, sin(A,—a)=1. 


Moreover, since there are only a finite number of roots of equation (2) lying in 
the interval 0 <A < X, where X, is the least value of 2, satisfying (12) and (13), 
or (12’) and (13), we may choose A, greater than 91/ and such that the inequality 
(11) is satistied for all the positive roots of (2) that are less than X,, provided we 
choose as the A, corresponding to each X* that root of (9) which lies nearest to 
r’. Hence we have, in view of (19), 


m 


(m=1,2,3,---), 


and the lemma is proved for the case of k >0. The proof for the case of k <0 
is quite analogous, and for k = 0 the proof is simpler. 

We will now show the relationship between equation (9) and equation (1) by 
proving the following lemma: 

Lemma 3. The positive roots of any equation of the form (1) are the same 
as the positive roots of a corresponding equation of the form (9). 


[> — 


1909] DEVELOPMENTS IN BESSEL FUNCTIONS 397 


We obtain easily from the asymptotic expansion of -/, (2) * 


2 2 1 1—4Y 2v+1 
(20) =a) [008 (a— + sin (x +6,(2)| 


where 
M, 
< xr (4>0), 


M, being a positive constant. 
We have for the derivative of J,() 


and on substituting this value in equation (1) we get 


(21) (vl + h)J,(%) =0. 


This equation may be thrown into another form by replacing J/,,,() by the 


v+ 


value obtained by changing v into v + 1 in (20) and also replacing -/,() by 


a value readily obtained from (20). We thus reduce (21) to the form 


2r 
-ige cos (x 4 


(22) 
+h)—1(4v+ 8v4+3). 
+ sin(a— 4 + 
where 
— | 
[e(A)|< 


M being a positive constant. 

Since the asymptotic expansions hold for all positive values of the argument, 
it follows that equation (22) is equivalent to equation (21) for all positive values 
of X, and hence the positive roots of any equation of the form (1) are the same 
as the positive roots of an equation of the form (9).{ 

We are now ready to prove the following theorem: 

THeorem I. Jf we represent by X,, the successive positive roots 
of equation (1), we have 


where C and K are positive constants. 


* Cf. LipscHiTz, Crelle’s Journal, vol. 56 (1859), pp. 193-196. 

+ Cf. GRay and MATHEWS, Treatise on Bessel Functions, p. 13. 

t There would seem to be an exception for the case/ 0. However, in that case (21) reduces 
by the substitution of (20) to equation (26) given below which is of the form of equation (9). 
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We prove this theorem first in the special case 7 = 0, in which case (1) 
reduces to 
(24) J,(r) =9. 


Sturm has shown * that for this case 


(25) lim —A,)=7. 
From this fact the first inequality (23) follows at once. 
In view of (20), equation (24) may be written for positive values of 


2 1 1— 4, 2 1 
(26) cos n)+ : sin(A— +6,(2)=0 52). 


Hence, however small the positive quantity 6, all the large roots of (24) lie in 
the intervals 


2v+1 
(27) T+ <6 (&=1, 2, ---). 


Since the distance between two successive roots of (24) approaches 7r as its limit, 
it is clear that one and only one of these roots lies in each of the more distant 
of the intervals (27). 


Precisely the same reasoning applies to the roots of the equation 


. Qv+1 
9 40S am» Ss => 
(28) cos (x 4 sin (a 4 0, 


since the distance between two successive roots here also approaches 7 as its 
limit, as we see from Lemma 1. 

Hence all the large roots of (26) coincide very nearly with the large roots of 
(28) and vice versa. 

Let us now denote by 2; that positive root of (28) which nearly coincides 
with the large positive root A, of (26).¢ Consequently A°’>C’n. By this 
inequality and Lemmas 1 and 2 we have for all large values of n 


| 


n n 


n 


For these large values of n we therefore have 


= | +> | + | A, — < a 


*Liouville’s Journal, vol. 1 (1836), pp. 174-175. 
{Notice that 7, is not necessarily the nth positive root of (28). However, 4, and 4,,, are 
successive positive roots of (28), and this is all that is essential. 
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If then we take K>2K'+ KH‘, we see that the second inequality (23) is 
established for all large values of x. It is clear that A’ may be increased so 
that it will also hold for the smaller positive values of n. Thus our theorem is 
proved when /= 0. 

Passing now to the case / + 0, we note first that by a theorem of SturM* 
between two successive roots of (24) there lies one and only one root of (1). 
From this, together with the fact that the first inequality (23) holds when 7=0, 
we infer at once the truth of this formula when 7 + 0. 

Reasoning now on equation (22) precisely as we reasoned above on (26), we 
see that however small 6 may be chosen, all the large roots of (1) lie in the 
intervals 


(29) 


<6 (k=1,2,---), 
and since these intervals for small values of 6 alternate with the intervals (27) 
we see from the theorem of Sturm last referred to that for large values of & 
one and only one root of (1) lies in each interval (29). 

That furthermore the large roots of 


2 2 
l(4v +8v+3) . (.- 


(30) cos (a in 


also all lie in the intervals (29) is obvious; that one and only one of them lies 
in each of these intervals follows from the fact that the distance between two 
successive roots approaches 7 as its limit, as we see from Lemma 1. 

Hence all the large roots of (30) coincide very nearly with the large roots 
of (22). 

From here on the reasoning is identical with the reasoning in the case / = 0; 
and our theorem is established. 

CoroLiary. Jf we represent by X,, r,, r,, the successive positive roots 
of equation (1), then there exists a constant q and a positive constant L such 
that 


L 
(31) |A,—ar—q|< (n ) 
Let us consider the series 


(32) 


n=1 


which by Theorem I is absolutely convergent. 
If we represent the value to which it converges by g, and its remainder after 
—1) terms by we have 


A, —(n—-1)74+ Rh, 


* Loe. cit., XIX. 
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and consequently 
=, K 

[A, 
“4 K n 2K 

< 5 dx = - =— 


n—1 


n (s=2)- 
Hence if we take g=9g,+.,—7, and L = 2K sufficiently large to satisfy 
the inequality (31) for n = 1, the truth of our corollary follows. 


§ 2. Summability and uniform summability of the development of an 
arbitrary function in terms of Bessel functions. 


We begin this section with a definition of summability and uniform summa- 
bility and the statement of a few general facts with regard to summable and 
uniformly summable series. 


Derinition. Jf u,(x), u,(x), are functions of x, and we let 
(33) 8. (x)= + 

(34) S(2)= = 


then if for a particular value of x, S,(x) approaches a limit U(x) as n 
becomes infinite, the series 
u,(x) + u, (2) 


is said to be summable to the value U(x); if S(x) approaches U(x) 
uniformly throughout a certain interval, the series is said to be uniformly 
summable.* 

It is a well known fact that every convergent series is summable, and every 
uniformly convergent series is uniformly summable, to the values to which they 
are convergent or uniformly convergent, respectively, but that the converse is not 
true. Moreover, the sum of any two summable or uniformly summable series is 
summable or uniformly summable, respectively, and its value is equal to the sum 
of the values of the original series. 

We will now proceed to establish some sufficient conditions that the series 


f af (x) 2) dex 
Fae 


where ,, A,, --- are the successive positive roots of equation (1), should be 


(35) 


 * The term summable is often used in other more general senses. With these we shall not be 
concerned in this paper. 


n=1 
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summable at every point of the interval 0 < x <1 at which the function f(x), 
which appears in the formule for the coefficients, is continuous or has a finite 
jump, and should be uniformly summable throughout every closed sub-interval 
lying within an interval of continuity of f(a). The series (35) is the ordinary 
formal development of an arbitrary function /() in terms of Bessel functions. 

We begin by proving two lemmas. 

Lemma 4. Jf we represent by X,,r,, --- the successive positive roots of 
equation (1), then the series 


(36) sin (A, 2 —a)dx* (p>4), 
n=1 


where a is any constant, will be absolutely convergent, provided $(x) is finite 
and integrable. 

Using the Law of the Mean and the corollary of Theorem I we obtain the 
formula 


1 1 
sin | (ur + q)u—a} dz 


L 
= $(x)(cos — nw — < = { \p(x)|\dx = 
where ZL, is a positive constant. We may therefore write 
Na 


nvo 


$(x) sin (A, — a) dx 
$(2) [sin — a) —sin + 2} ]de 
l | 


a $(x) sin {(nr + q)x— a} 
But by Theorem I 


(37) 


~ Cen'te? 


and hence the series whose general term is the first term on the right hand side 
of (37) is convergent. If now we can show that the series which has for its 
general term the second term on the right hand side of (37) is convergent, it 
will follow that the series (36) is absolutely convergent. 


*We have taken the sine as the trigonometric function appearing in the integrand. That the 
lemma holds equally well for the cosine is seen by making the change of variable 


3° 
~ 
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We have 


1 
(38) sin + — a} do =u, +0, 


where 
l 1 
(x) cos (qx — a) sin nradx, = (x) sin (gx —a) cos nrradz. 
2/0 0 


These quantities «, and v, are constant multiples of Fourier’s constants belong- 
ing to the functions cos (qu —a) and $(x) sin (gv — a), respectively. 
Therefore, since these functions are finite and integrable, it follows from a the- 
orem due to DE LA VALLEE Poussin that 


x 
9 
n=1 a= 
are convergent.* Hence the series 


w 


(39) 


> 
V 


are convergent. +} 
But we have from (38) and Theorem | 


|e, | lo. | | 


1 
(40) f = he = Ce ne’ 
Since the two series in (39) are convergent, it follows from (40) that the series 
whose general term is the left hand side of (40) is convergent, and our lemma is 
proved. 
Lemma 5. If the general term of the series (35) is written in the form 
9 
cos (qx — a) cos nx V cos (qa — x) cos dx 


(41) 


vl 
+ ——sin (qx — a) sin V.af (x) sin( — a) sin + v, (x), 
V x eo 


then, provided the function f(x) is finite and integrable, the series 


(2) 


n=1 


will be uniformly summable in every interval 
(42) c=r=1 

*Annales de la Société Scientifique de Bruxelles, vol. 17 (1892-3); also 
BocHER, Annals of Mathematics, ser. 2, vol. 7 (1906), p. 107. 

t BOcHER, loc. cit., p. 108, lemma. This lemma is there proved only for the case p —1 but 
the proof applies equally well to the general case p > }. 


Iv | 
n 
| 
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for which 
and if we write 


= xv (x) (0<zr=1), 
6,(9)=lim (x) 


the series 


Lz, (2) 


n=1 


will be uniformly summable and its terms will be continuous in the intervat 
(43) 0=27=1. 


It is obvious that it will be sufficient to prove only the second part of the 
above lemma, since the first part is an immediate consequence of the second. 
We will show then that the facts stated hold with regard to the series whose 
general term is 0, (2), and it will follow at once that the facts stated hold also 
with regard to the series whose general term is v,(#). 

In view of the asymptotic expansion of J,() [see equation (20)] we may 
write 

A ¥ (A) 
(44) JS(A)= a (A — a) + 


where A and a are constants and we use (2d) to represent any function of X 
that is continuous for all positive values of X, remains finite for all such 
values of X%, and approaches a finite limit when X approaches zero. It is 
obvious from this definition that in different formulas, and even in different 
parts of the same formula, (2) may have a different meaning. 

If we set Aw = y we obtain 


But we have * 


and hence we obtain from (44) and (45) as a value for the denominator of the 


fraction in (35) f 


(46) [2 = + ) 


where y,() is used to represent any function of X that remains finite when 


* Cf. BYERLY, An Elementary Treatise on Fourier’s Series, etc., p. 224. 
+ For the sake of simplicity we use 7 instead of 7.,. 


r=+0 
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X increases indefinitely. It is apparent from this definition that a (2) is 
always a but that a may not be a 

Turn now to the numerator of the fraction in (35). If we use ¢ to denote 
any constant lying in the interval 0 <#<1, and set Ax = y, we obtain for 
A>c 


a1) f "af (de) de = 


Moreover, we have from the asymptotic expansion [see equation (20) ] 


(48) = | cos (ae — a) sin (Aw — a) + |. 


From (48) we obtain 


A A 


(re: 
V Yr 


a A’ 
= Vax f(x) cos (Ax — + xi 


(49) 
art (Aw)! 


V 


Combining (47) and (49) we get 
el 4 1 
xf (x) S,( Ax) dx = Vax f(x) cos (Ax — a)dx 
0 V 0 
(50) 
7 sin (Ae — + ) 


Hence, for the general term of (35) we obtain by combining (44), (46), and 
(50) 


1,008) f 
(51) 


A 
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where 
el 


af 
Z 


f(x) é 
sin (Ax — 1; 


If we subtract from the right hand side of (51) the two terms 


9 1 
(52) cos — a) f Va f(x) cos (Ax — a) dx, 
0 


2A’ 


(53) cos (Aw — a) sin (Ax — a) dx, 


AXV 2 


the quantity we have left may be written in the form 


4) A?’ ( ) 


where W(a, X) is used to represent any function of « and X which for all fixed 
values of d is continuous in the interval (48) and which remains finite when r 
increases indefinitely and x varies in the interval (43). 

We have now broken up the general term of the series (35) into the sum of 
the three expressions (52), (53), and (54). It is obvious that (53) and (54), 
when each is multiplied by and properly defined* for «= 0, are con- 
tinuous in the interval (43); we shall show further that each of these expres- 
sions when multiplied by x? forms the general term of a series which is uniformly 
summable in that interval. Then we shall prove that if we subtract from (52) 
a number of expressions which when multiplied by «3 and properly defined for 
a = 0 are continuous in the interval (43) and form the general terms of series 
that are uniformly summable in that interval, the expression we have left will be 
equal to the sum of the first and second terms of (41). Our lemma will thus be 
established. 

We will begin with the proof that (54) when multiplied by «! forms the gen- 
eral term of a series that is uniformly summable in the interval (43). For val- 
ues of A so great that 

Wi(A)) 


| 


where the y,() occurring in this inequality is the same as the w,(2) in the 
denominator of (54), the expression (54) when multiplied by ? is less in absolute 


value than 


*T. e., defined to have the value it approaches when x approaches zero. 
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+ 4 | 


Since x(a) is finite and integrable in the interval 0 = # = 1 it follows from 
Lemma 4 and the first inequality in (23) that the series whose general term is 
(55) is absolutely and uniformly convergent and therefore uniformly summable 
in the interval (43). 

We will now show that (53) when multiplied by 2! is the general term of a 
series that is uniformly summable in the interval (43). If we set M=2A'/A, 
(53) when multiplied by x! becomes 


Mz cos (Ax — 
(56) Mx cos(Axwn— a) (*' f(a). 
¢ sin (Ax a)dx. 


The second term in (56) is equal to 


Mz cos — a) 
| F(a) sin (Ax — a) dx, 
where we define 
7 ( ) <= <= \ 7 <= 
F(2)= (c=z=1), F(x)=90 (0 S2<ce), 

and hence it follows from Lemma 4 that it is the general term of a series that 
is absolutely and uniformly convergent and therefore uniformly summable in the 
interval (43). 

We have next to examine the series whose general term is the first term in 
(56). It follows from the corollary of Theorem I that 


and consequently 
x | 
(58) cos (A, 2 — a2) = cos + — a | cos n) 


—sin [(nr+q)x—a] sin 


(A,) 
n 


n 


= cos 


In a similar manner we may also obtain from (57) 
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e, A 
By the use of (58) and (59) and the first inequality in (23), the first term in (56) 


reduces to 


COS q) a | nor + — a] de + 
V2 


n 


(60) 


n 


and since the second term in (60) is the general term of a series that is abso- 
lutely and uniformly convergent and therefore uniformly summable in the inter- 
val (43), we have only to examine the series whose general term is the first term 
in (60). 

Using the value of X, given by (57), we see that if we add to the first term in 
(60) the term 


+ Mz cos + q)x— a] 


na (um + 


which is easily seen to be the general term of a series that is absolutely and 


sin + q)x —a]dz, 
o 


uniformly convergent, and therefore uniformly summable in the interval (43), 
we obtain the expression 


Mx cos f(x) 


61 i lx. 
(61) sin [ + q)x— a] da 
If now we can show that the series which has (61) as its general term is uni- 
formly summable in the interval (43), we shall have proved that the series which 
has (53) multiplied by x? as its general term is uniformly summable in that 
same interval. 

By a familiar trigonometric formula, (61) reduces to the form 
Me f(2#) 
J, x’ 


{sin + q)(x +2’) — 2a] de’; 


which, if we set K = M/2z7, may be written in the form 
sin +2 ) te 


ke [22 cos [g(a + a’) — 2a] - 


+ Ke [7 sin [q(a + — 22 


($3) 
Ke eos g(x — x’) — dx 

0 V2’ n 

n 


Trans. Am. Math. Soc. 27 


A, ) 
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We are going to show that each term in (62) is the general term of a series 
that is uniformly summable in the interval (43). 
If we represent by s (z) the sum of the first n terms of the series 


=. sin nz 
(68) 

n=1 
we have for values of z such that sin (z/2) + 0* 


8, (2) 


(Qn4+1lyz 1 (2v + 2p 1)z 


2 sin 


In an analogous manner we obtain for the series 


=. COS NZ 
(65) 


n=1 
the relation 


8,(z) — §,(2) 


1 Al 1 1 1 
(67) tp 1 
<- (p=—n+1, n+2, ---). 


2 sins - (n+1)- sin 5 


In an analogous manner we obtain for the series (65) by use of the relation (66) 
the same inequality (67). 

If now we represent by s (2) the sum of the first n terms of the series whose 
general term is the first term in (62), we obtain by means of (67), for all values 
of x in the interval 0 << x = 1, the inequality 


(x 2) —3,(2)| Ke f (# +2 2a] di! 


v=n+l 


*Cf. page 110 of BOcCHER’s article referred to above. 


2 sin 5 
In view of (64) we have for the series (63) 
| 8,(2) — 8, (2) | 
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= x+a’ V2 ‘ ~n+1 

(p=n+1,n+ 2, ) 


where K;, is a positive constant. For x = 0 the last inequality obviously holds, 
since each term of the series we are discussing is equal to zero. Hence it fol- 
lows that the series whose general term is the first term in (62) is uniformly 
convergent and hence uniformly summable in the interval (43). In a precisely 
analogous manner we may establish by use of the relation (67) the fact that the 
series whose general term is the second term in (62) is uniformly summable in 
the interval (43). 

We will consider next the series whose general term is the fourth term in 
(62). If we represent by s(x) the sum of the first n terms of that series we 
have, in view of (67), 


(2) —8,(#)|S Ke dx’ 
0 V v=a+l1 n 
= Kx ( x’) | | sin q (x ) | ; dx’ A, 
a n+1 


(xn +1)-|sin 
(p=n+1,n+2, 


for all values of x in the interval (43), A, being a positive constant. Hence the 
series whose general term is the fourth term of (62) is uniformly convergent and 
therefore uniformly summable in the interval (43). 

We have finally to consider the series whose general term is the third term 
in (62). This term may be written in the form 


Kx cos gz ) cos gar’ sin (x2 — ) de’ 
(68) 


sin nw — 2’ 
— ( id 


nr 


+ Kz sin gx | f(x ) sin 
By means of an integration by parts we may reduce the first term of (68) to the 
form 


Kx cos ga — ) cos ga’ dx 
c 
(69) 
+ 7Kx cos qx f 2’) cos — x’) dz’, 
c 


where we have set 


(70) F(z’) = cos (eZ 220). 
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If we represent by s(x) the sum of the first n terms of the series whose gen- 
eral term is the first term in (69), we have, by means of (67), for all values of 
x in the interval O<v=1, 


K, 
re ~n+1 
2 


cos ger [ cos ga'da’ 
(nm +1) -sin 


(p=n+1,n+2,---), 


where , is a positive constant. For « = 0 the last inequality obviously holds, 
since each term of the series that we are discussing is equal to zero. Hence it 
follows that the series whose general term is the first term in (69) is uniformly 
convergent and therefore uniformly summable in the interval (43). 

Let us now consider the series whose general term is the second term in (69). 


If we set 
(71) u, (2%) = 008 (2 — x')dx’, 
where i 
x(x’) = F(x’) (0S2’Sce), 
(72) 
x(#)=90 (e<a/S2), 


the second term of (69) may be written in the form 

(73) — cos 

However, it follows from (71) that u,(a) is the general term of the Fourier 
development corresponding to the function x(a), and since we know from (72) 


and (70) that y(<’) is finite and integrable in the interval 0 =a2=2, it follows 
from FEsé£R’s work * that for the series whose general term is u,(«), 


(74) |S (s-=1, 3, -*-; @Se359), 
where S (a) is defined as in formula (34) and J is a positive constant. If 


now € is given, positive and arbitrarily small, we have from (74) for the series 
whose general term is (73) or the second term in (69), 


(75) — S,(x)| S,(2)| < cos = 


where we have set B = 27K N. Moreover it follows from another theorem of 
Fesér’s} that the series whose general term is u,(x#), and hence the series 
whose general term is (73) or the second term in (69), is uniformly summable in 


the interval 


* Loe. ecit., p. 60. 
t Loe. cit., p. 60. 


€ 
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since x(x) is continuous in an interval including this interval. Hence we may 
choose m so large that we have for the latter series, 


(76) | — <e (p>nZm; 


Combining (75) and (76), we see that the series whose general term is the second 
term in (69) is uniformly summable in the interval (43). 

Since both the terms in (69) are the general terms of series that are uniformly 
summable in the interval (43), it follows that the series whose general term is 
the first term of (68) is uniformly summable in that interval. In an entirely 
analogous manner it may be shown that the series whose general term is the 
second term in (68) is uniformly summable in the interval (43), and hence it 
follows that the series whose general term is (68), or the third term in (62), is 
uniformly summable in that interval. 

We have now shown that each term in (62) is the general term of a series 
that is uniformly summable in the interval (43). Hence (62), or (61), is the 
general term of a series that is uniformly summable in that interval and conse- 
quently, as we have pointed out before, (53) when multiplied by ? is the gen- 
eral term of a uniformly summable series in the same interval. 

It remains now to consider the expression (52). It follows from (58) that if 
we subtract from (52) a term of the form 


(x, r,) 


( ) NV 


1 
Vaf(x) cos — a)dx, 
0 
which when multiplied by 2? and properly defined for x = 0 is obviously continu- 
ous in the interval (43), and by Lemma 4 is the general term of an absolutely 
and uniformly convergent and therefore uniformly summable series in that same 
interval, we have left the expression 


1 
(78) eos + — a] Vaf(#) cos (A, 2 — a)dx. 
0 


Moreover we obtain from (58) 


,) 


sin 


cos (A, = cos + a] 


and hence (78) reduces to 


1 
7x 008 + q)x—a] Vaf(x) cos + q)x—a} dz 
Vz 0 
(79) 


nV 


¥ 


WV 


1 
cos [(n7+q)x—a] f (2) sin dx+ 
0 


The last term in (79) when multiplied by x! and properly defined for x = 
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is obviously continuous in the interval (43) and is the general term of a series 
that is uniformly convergent and therefore uniformly summable in that interval. 
The second term in (79) when multiplied by x! and properly defined for x = 0 
is obviously continuous in the interval (43); that it is also the general term of a 
uniformly convergent and therefore uniformly summable series in that interval 
may be seen by referring to the proof of Lemma 4. It is apparent from the 
latter part of that proof that a series whose general term has the form 


(80) $(x)sin [(ur +q)x—a]dzx, 


where $(«) is finite and integrable, is absolutely convergent. The second term 
of (79) when multiplied by «x! consists of a factor of the form (80) multiplied by 
a factor that remains finite for all values of X, and all values of x in the interval 
(43), and hence this term when multiplied by x! is the general term of a series 
that is absolutely and uniformly convergent and therefore uniformly summable 
in the interval (43). 

It only remains to consider the first term in (79) which may be broken up 
into the following sum of four terms: 


9 el 
(81) cos (qx — a) cos nore | V xf (x) cos (qu — a) cos nx dx 
Vax 
9 1 
— — cos (gx — 2) cos V xf (x) sin (gx — @) sin nx dx 
Va 


0 


sin (qx — a) sin V xf (x) cos (qx — a) cos dx 
V & 
1 
+ —=sin (qv — a) sin ner V (x) sin (gx — 2) sin 
Vz 


The first and fourth terms of (81) are identical with the first and second terms 
of (41). The second and third terms when multiplied by x? and properly defined 
for « = 0 are obviously continuous in the interval (43); if we can show further 
that they form the general term of a series that is uniformly summable in that 
interval, our lemma will be completely established. 

If we use a’ as the variable of integration, the second and third terms of (81) 
when multiplied by 2? may be written in the form 

1 
of sin [ + x) — 2a] sin + 2’) de’ 

(82) 


—x} g(x — 2’) sin —2’)dz’. 


1909] DEVELOPMENTS IN BESSEL FUNCTIONS 413 


Since, however, 
sinz — sin(n +1)z + sin nz 
sin nz = —— - 


2 — 2 cosz 


we obtain for the series whose general term is the expression (82) 


( x) 


1) sin 2’)—sin ( 
(83) 4n sin? 7 


9 


sin — sin (2 4 1 
0 


9 


4n sin? 7 


where S (x) is defined as in formula (34). 

We are going to show that as m becomes infinite each of the terms on the 
right hand side of (83) approaches a limit uniformly for all values of # in the 
interval (43) and consequently that the series whose general term is (82), or the 
second and third terms in (81) multiplied by x, is uniformly summable in that 


interval. 
Consider first the second term on the right hand side of (83). It may be 
written 
: nt sin Vai f(z’) sin — ©) ax! 
sin sin T 
(84) 


4 — (n +1) Vx—27 V2—2' 


sin — sin — 


x sing(a—a’) sin dx 
0 


The integrand of the integral in the first term of (84) consists of a function of 
x and 2’ that is continuous in both variables multiplied into a function of x’ that 
is finite and integrable in the interval of integration. Hence that integral 
defines a function of x that is continuous and therefore finite in the interval 
(43), and consequently the first term in (84) approaches a limit uniformly for 
all values of x in the interval (43) as x becomes infinite. 
The factor 
(85) sin (n + 2") 


(m+ 


which occurs in the integrand of the second term in (84) approaches zero uni- 


n=1 
1) sinz —sin(n+1)z 
n 2n(1— cosz) 
| 
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formly for all values of 2 and x’ as n becomes infinite. For, given € positive 
and arbitrarily small, we have for values of x and 2’ such that 


é 
2? 


- 


(n+1)V 
Moreover, we can choose m so large that for all values of x and 2’ such that 
e 
|= 


we have 
sin (n + 1)r(x—~2’) 


(n+1)Va—2 


(n=m). 


~ 


We see then that the integral in the second term in (84) is in absolute value 


not greater than 


 sin(n+1)r(2—2’) dx’ 


(86) Jo sin T sin 
el dx’ 
< Ke — (n2m; 05221), 


where KX is a positive constant. But 


(87) f dat on <2f da 
o |) la| 


| | J | 


Combining (87) with (86) we see that the integral in the second term in (84) 
approaches zero uniformly for all values of « in the interval (43) as x becomes 
infinite, and hence the second term itself approaches zero uniformly in the same 
interval. 

The first term on the right hand side of (83) may be broken up into the sum 
of the following three integrals : 

x ‘sin [g(x +2’) — 2a] 


; 
Va' f(x’) ( 


‘sin [g(a + 2’) — 2a sin +a 
(88) ; J V f(x’) - ( ) de’, 


2 
sin 7 - n sin 7 


. 2 2 
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provided that we can show that each of these integrals converges for all values 
of «x in the interval (43). For this purpose we must consider the values of q. 
By comparing equations (20), (44), and (22) we see that (1) may be written 

—1(47 + 8v 43) 


(A —a)— Alsin (A — a) 


Consequently, except in the case / = 0, the large roots of this equation coincide 
more and more nearly with the values 


a+tkr, 


where & takes on large positive integral values. In the exceptional case / = 0 
they coincide more and more nearly with 


g t+ athr. 


Since, by (31), we see that 2, coincides more and more nearly with g + nz, it 
follows that according as the / of equation (1) is or is not zero, 


T 
or q=at+hk,r, 


where k, and k, are positive or negative integers or zero. Consequently the 
expression 
sin [g(x +2’) — 2a] 
(89) 
sin 

approaches a finite limit as 2 + «’ approaches the limit 2, and hence the inte- 
grals in (88) converge for all values of x and all values of x in the interval (43). 

It now remains to be shown that the expression (88) approaches a limit 
uniformly for all values of x in the interval (43), as m becomes infinite. The 
first term of that expression is independent of x and hence need not be consid- 
ered further. The integral in the second term has for its integrand the product 
of a function of 2 and 2’ that is continuous in both variables and a function 
of x’ that is finite and integrable. It therefore defines a continuous function of 
#, and consequently the second term approaches zero uniformly for all values of 
x in the interval (43), as n becomes infinite. 

The same is true of the third term of (88). For, since the expression 


sin(n +1)7(x+~2’) 


e+ 
nsin 


remains finite for all positive integral values of n and all values of x and 2’, it 
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follows that if € is given positive and arbitrarily small, we can choose 7 so small 
that 


sin 7 n sin 
2 2 


(a=1, 8, 3, 
We ean then choose 6 so small that 


(91) if "sin [ g(x +2) 2a | Vai f(x’) +1)r(a + a ) 


sin n sin 


Finally, 7 and 6 being fixed, we can choose m so large that 


(92) nf e+e Va f(x) da’ |< 
0 7 - sin — 


(n=m; 


= 


Combining (92), (91), and (90) we see that the third term of (88) approaches 
zero uniformly for all values of x in the interval (43) as n becomes infinite. 

We have now shown that as n becomes infinite each of the terms on the right 
hand side of (83) approaches a limit uniformly for all values of x in the interval 
(43) and therefore it follows that the series whose general term is (82), or the 
second and third terms of (81) multiplied by x’, is uniformly summable in the 
interval (43). Hence, as we have pointed out above, our lemma is proved. 

We are now prepared to prove the following theorem: 

TueoreM II. /f in the interval 0 =x =1 the function f(x) is finite and 
integrable, then the series (35) will be summable at every point of the interval 
0<x<1 at which f(x) is continuous or has a finite jump; and will be 
uniformly summable in every closed interval lying in that interval which does 
not include a point of discontinuity of f(x). 

Since f(a) is finite and integrable, 


Vax f(x) cos (qu — a) and V « f(x) sin (qx — a) 


are finite and integrable, and hence it follows from a theorem due to Fesér* 
that the two series whose general terms are the first and second terms of (41), 
respectively, are summable at every point of the interval 0 < x < 1 at which the 
function /(#) is continuous or has a finite jump, these series being essentially 
equivalent to the cosine development and sine development, respectively, of 
Vax f(x) cos —a) and Vx f(x) sin(qgz—a). We know from Lemma 5 


* Loc. cit., p. 53. 


1). 
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that the series which has for its general term the third term of (41) is uniformly 
summable throughout the interval (42), and hence, since the c of that interval is 
arbitrary, is summable at every point of the interval 0<2<1. Consequently 
the series (35), whose general term is the sum of the three terms occurring in 
(41), is summable at every point of the interval 0 <2 <1 at which f(x) is 
continuous or has a finite jump. 

We know further from Fres&ér’s work * that the series having for their general 
terms the first two terms in (41) are uniformly summable in every closed inter- 
val which lies in the interval 0 < a < 1 and does not include a point of discon- 
tinuity of f(a), and it follows from Lemma 5 that the series having for its 
general term the third term in (41) is uniformly summable in any such interval. 
The series (35) is therefore uniformly summable in any such interval and our 
theorem is completely established. 

Coro.tiary. If we define S (x) as in formula (84), we have for the series (35) 


|S(2)|<M  (n=1,2,3,---; 0<eSx=1) 


where M is a positive constant. 
For each of the two series whose general terms are the first and second terms 
of (41) we have from FEs£R’s work { 


(93) <M, 


where JV, is a positive constant; for the series whose general term is the third 
term of (41) we have from Lemma 5 


(94) | M, (n=1, 


where JV, is a positive constant. If we combine (93) and (94) the proof of our 
corollary follows at once. 

Before proceeding to the proof of the next theorem, we wish to make the fol- 
lowing convention. Let 


(95) $(x) = cos’ (gx — a) f(x) + sin’ (qx — a) f(x) + yv,(2), 


where v, has the same meaning as in (41), and the series is to be evaluated as a 
summable series. The function ¢(a) we will designate as the function associ- 
ated with the series (35). It will be seen that ¢(a) is thus defined at all 
points where f(a) is defined, whereas the series (35) is not necessarily sum- 
mable at points where f(a) is discontinuous. At points where /(2) is contin- 
uous, however, it is readily seen by a reference to Lemma 5 that the series (35) 
is summable to precisely the value ¢(2). 


* Loe. cit., p. 60. 
T Loc. cit., p. 60. 
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THeoreM III. Jf we multiply the series (35) by any function x(x) that is 
finite and integrable in the interval 0 = x = 1, and integrate the resulting ser- 
ies term by term from &* to 1, we shall obtain a series which is summable to 


the value 
(96) p(x)x(x)dx 
Js 


where $(x) is the function associated with the series (35). 

Since the series whose general terms are the first and second terms of (41) 
differ from Fourier developments of finite and integrable functions only by fac- 
tors which are continuous functions of « in the interval 6 = 2 = 1, we know from 
a theorem due to DE LA VALL&E Poussin ¢ that we can multiply each of them 
by any finite and integrable function (a) and integrate the resulting series 
term by term from 6 to 1, and that the series thus obtained will converge to 
the value of the integral of the product of y(2x) and the function of which the 
series is a development multiplied by the function which appears as a factor of 
the series. 

Since, from Lemma 5, the series whose general term is the third term of (41) 
is uniformly summable in the interval (42), we can multiply it by any finite and 
integrable function (a) and integrate it term by term from 6 to 1, and the 
resulting series will be summable to the value of the integral of the product of 
x(«) and the function to which the series is summable. Thus our theorem is 


proved. 
$3. Value of the development. 


The series (35) is usually given as the formal development of an arbitrary 
function in terms of Bessel functions, but there are cases in which this series, 
even if it is convergent or summable, will not have the desired value, an extra 
term being necessary in order that the development may represent the function 
developed.{ However, in these cases the extra term appears naturally as part 
of the formal development if we keep iy mind the method by which the formal 
development is commonly deduced. 

In the problems of Mathematical Physics in which we wish to develop an 
arbitrary function in terms of certain normal functions, the latter functions 
appear as solutions of some ordinary differential equation. In the case in which 
the normal functions are Bessel functions the differential equation is usually 


(97) 


dx? x dx 
* The quantity 4 is any positive quantity less than 1. 
+Cf. Annales de la Société Scientifique de Bruxelles, vol. 17 (1892-3). Cf. 
also BOCHER’s article referred to above, p. 118, footnote. 
{Cf. Dini: Serie, di Fourier, and BRIDGMAN, Philosophical Magazine, vol. 16 (1908), 
p. 947. 
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where v is a constant = 0 and X is a real constant (which may obviously be 
taken = 0) to be determined in such a way that (97) has a solution, not iden- 
tically zero, which remains finite when « = 0, and at the point x = 1 satisfies 
the condition 


(98) [ Zu’ (a) + hu(a)],_,=90 


where / and / are givea constants not both zero. 
If in (97) A + 0, the most general solution which remains finite at x = 0 is 


(99) u= CJ,(rAx), 
and in order that this should satisfy (98) it is necessary and sufficient that 
(100) + = 9. 


This is the equation for determining the values of X greater than zero which we 
have to use. When, however, \ = 0, formula (99) no longer gives the general 
solution of (97) which remains finite when « = 0, this solution then being 


(101) u = Cr’. 
The latter solution satisfies (98) when and only when 
(102) vw+h=0. 


It is, therefore, only in this case that the value X = 0 needs to be considered 
at all. In all other cases we develop our function in the form 


f(x) = DAA 


where A,, A,, --- are the successive positive roots of (100), and the coefficients 
A are to be determined by the ordinary formulas. 

In the exceptional case (102) the A’s greater than zero are still determined 
by (100), which now reduces to 


(108) —1J,,,(%) =0.* 


But we also have to use the value } = 0, and thus we are led to the develop- 
ment 


(104) f(2) = + SAS, (2,2), 


n=2 


where A,, are the positive roots of (103). Here not only are VxJ, (A, x) 
and VxJ,(A,x) orthogonal to each other, i. e., 


1 


~ *Cf. for instance BYERLY, loc. cit., p. 223, formula (6). 
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but x’*?is also orthogonal to each of the functions VxJ,(A, x), i. e., 
vl 
i] J (vA 2) da = 0. 


Accordingly here also the coefficients can be determined in the ordinary way. 
For instance, to get A, we multiply (104) through by «’*' and integrate term 
by term from 0 to 1, thus obtaining * 


A 
xv’ A, = +. 9? 


l 
A, =(2v+ 2) 


If then, in developing our arbitrary function, we use a// of the normal fune- 
tions that satisfy (97) and (98) we shall naturally have an extra term in the 
case in which (102) is satisfied. Hence, what appears to be an exceptional case 
is really not exceptional at all, save in regard to the notation ordinarily used 
for expressing the solutions of (97) in terms of Bessel functions. Moreover, it 


so that 


is possible to adopt a notation which avoids the awkwardness of making a special 
case of the solution (101). For, if we set 


F,(4,2)=(5) T(y +i +1) 
so that 
(105a) J, (Av) = (A, 
then 

u= CF,(r, x) 

will be the most general solution of (97) which remains finite when a = 0, not 
only when A> 0 but also when A = 0, since in this last case it reduces to the 
form (101). 


We now consider the series 


(106) 2 —=) A 


l 
n=1 f x) ]?dx n=1 


where A,, A,, A,, --- are the roots, positive or zero, of the equation 


(107) Fe) + | = 04 


* It is understood of course that we are using here merely the formal device ordinarily employed 
for finding the coefficients of the development without entering into the question of whether or 
not we have a right to integrate the series term by term. 

Tt The positive roots of this equation are obviously identical with the positive roots of equation 
(1), and hence all that we have proved for the latter roots holds for the former ones. 


4 
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arranged in increasing order of magnitude. The equation (107) will have a 
zero root when and only when (102) is fulfilled, and thus the series (106) will 
give the formal development of the arbitrary function f(2) for all cases. This 
series will be identical with (35) except when (102) is fulfilled, and in that case 
will differ from it only in having an extra term. Since, however, this term will 
be finite and continuous ic will not affect the summability or uniform summa- 
bility of (85) and all that we have proved in the previous section with regard 
to that series will hold in all essential respects with regard to the series (106). 
Thus if we write the general term of (106) in the form 


1 
cos (qx — a) cos nora Vx f(x) cos (gu — 2) cos 
0 


V 
(108) 9 
sin sin nore V (x) sin (gv—a)sinnrxdx+w, (x), 
0 


the series 


n=1 


will be uniformly summable in the interval (42), and if we define the function 


(109) (a) = cos’? (qr—a) f(x) + sin? (qa —a) f(x) + 


to be the function associated with the series (106), Theorem III holds without 
further change if we replace the series (35) by the series (106) in the statement 
of that theorem. Moreover ¢(x) will have the value to which the series (106) 
is summable at every point at which f(x) is continuous. 

We now wish to show that the series (106) is summable to the value 


(110) +0) —0)] 


at every point of the interval 0 <a <1 at which /(2) is continuous or at which 
it has a finite jump, provided the latter points lie within an interval at every 
other point of which f(x) is continuous. We will begin by proving some 
lemmas. 

Lemma 6. If f(x) is finite and integrable and $(x) is the function asso- 
ciated with the series (106), we have 


1 1 


where X,, Az, As, +--+ are the roots, positive or zero, of equation (107), arranged 
in increasing order of magnitude, and F(X, x) is defined by equation (105). 


i 


422 C. N. MOORE: SUMMABILITY OF [October 


We have, from the modification of Theorem III just referred to, the follow- 
ing equality * 


1 
xp(x)F(r,, «)dx 
(112) “° 1 
Js 3 


Moreover F',(X,, x) and F,(A,, satisfy differential equations of the form 
(97). From these we readily deduce, if we use accents to denote differentiation 
with regard to x, 
(2 — )aF,(A,, (A,, 


n? 


d 
= dz” [F(A 


d 9 7 9 
and consequently when X, + A, 


vl 
J (r,, 7) F,(A,, «) dx 
(113) 


== 

It follows without difficulty from equations (105a) and (44) that at the upper 
limit the quantity in square brackets on the right hand side of (113) may be 


written 
( ) 
At the lower limit it may be written - 
+4 
if 8 lies in the interval 
c 
0 < 8 ~, ¥ 


< It is understood of course that the series on the right hand side of (112) is to be evaluated 
as a summable series, as we have yet no reason to suppose that it is convergent. 
t This last reduction follows from the formula 


Fi(4,2)=—#P x) +-F,(4,2). 
Cf. BYERLY, loc. cit., p. 223, formula (6). 


n 
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where c is any positive constant < 1, since we have from (105) 


If 5 > c/X, it may be written 
¥,(A,) 


since for such values of 6 we have from (105a) and the asymptotic values of 


J, (x) and J,,, (x) 


9) Vx, arr (A, BL Vor 

¥,(A,) 

v+l1 ( n ) V +3 


Finally if 6 = 0 the expression in square brackets on the right hand side of 
(113) vanishes at the lower limit. Hence we have 


M 
(114) f (r,, 2) F,(A,, dx (nk, OS5<1), 
where 1 is a positive constant. Furthermore it follows from (46), (50), and 
(105a) that 
Zarti 
(115) cos — a) dx + (2, ). 


Consequently, in view of (114) and (115), we have for the general term of the 
series on the right hand side of (112) 


1 
4, F(A,; x )dx 
8 


2Mr, 


The second term on the right hand side of the above inequality is, in view of 
the first inequality in (23), the general term of a convergent series, and it follows 
from Lemma 4 that the same is true for the first term. Hence the series on 
the right hand side of (112) converges uniformly for all values of 6 in the 
interval 0 = 6 < 1, and moreover it converges to the value of the integral on 
the left hand side for all values of 6 > 0, since it is summable to that value. 
Since the convergence is uniform, we may let the 6 of equation (112) approach 
zero. But in view of the fact that for the different values of X, that enter, the 
functions V. a F, ( X,, 2) are orthogonal to each other, all the terms but one on 
the right hand side of (112) drop out, and our lemma is established. 


Trans. Am. Math. Soc. 28 
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Lemma 7. The function $(x) associated with the series (106) does not 
become infinite to an order higher than x~i as we approach the point x«=90. 

By definition the function ¢(a) is equal to the sum of two functions that are 
finite and integrable in the interval (43), plus the series 


(116) »,(#)- 
n=l 
By Lemma 5 the series (116) is equal to the product of x~ into the series 


v, (x). 

But by the same lemma this latter series is uniformly summable and its terms 
are continuous in the interval (43). Hence it defines a continuous function in 
this interval, and consequently the series (116) defines a function that does not 
become infinite to an order higher than x~? when we approach the point x = 0. 

Our lemma is therefore proved. 

TueoreM IV. Jf we have a function x(x) which in the interval 0 Sx =1 
is finite everywhere save in the neighborhood of the point x =9, and if the 


1 
( x ) dx 


converges absolutely, and if furthermore 


integral 


(117) = 0 (n=1, 2, 3,°--) 


where X,, r,, rs, «++, are the roots, positive or zero, of equation (10T) arranged 
in increasing order of magnitude, then y%(x) is zero at every point at which it 
is continuous. 

We know from some theorems due to Din1* that we can develop each of the 
functions 


(118) sin kore (k=1, 2, 3, ---) 
in a series of the form 
(119) 2%) + a,F,(A,, + a, + 


and that the series will converge to the value of the function developed at every 
point of the interval 0 <2 <1. Weare going to show further that for the 
functions we are considering the developments will converge uniformly through- 
out the interval0 =2=1. 

* Serie di Fourier, vol. 1, pp. 246-269. DINI treats a very general class of functions of which 
the functions (118) are a special case. 


t There is an error in Drn1’s formule for the special case in which 2, = 0; he omits the factor 
F, (4,, 2) from the first term of the series. 


| 
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We will consider first the coefficient of the general term of the series (119) 
and for simplicity write X for ,.. We have then 


1 
a sin dx 
(120) 


n 


1 
F(A, x) ] "dx 
It follows from equation (46) and equation (105) that 


l A? x 


We may therefore choose d so large that 


(121) ja,|<- 


But we have 


1 
1 
(122) 2? sin 2) de = + 
0 0 
A 


and if we set Ax = y 


1 
sin kraF’,(r, «)dx. 


1 1 


1 A 
“x? sin «)dx = sin kore J, (Ax) dx 
(123) 


1 
y sm x J,(y) dy. 
Moreover we have from equations (49) and (105a) 


1 
x sin krraF’,(r, x) dx 


rr. 
= sin cos (Aw —a)dx+ f xt sin sin (Ax —a)dx+ ) 
*J 1 
(124) A A 
A 1 A rel 
“Jd 
if x! sin cos (Aw—a)dax— if sin sin (Aw—a)da+ ) 


But we have from equation (57) for ’ = A, * 


cos (Ax — a) = cos (ume + qx — a+ 


*In the special case in which 2, = 0 we must take 2 == 7, , for (125) to hold true. How- 
ever, since we are dealing here only with a question of convergence, this is unessential. 


| | 
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ay (A (r 
(125) 
x 
. 
— sin + ga — a) sin 
Moreover 
el 
| al sin kare: cos gx—a) dx x4 sin cos (gr—2) cos dz 
= wv! sin sin (gx — a) sin dx ) 


and from an integration by parts 


ay, 
x? sin sin + ga — a) sin vil ») d. 


= |- x! sin kara cos + gu — a) sin ( | 
+ q 0 
(127) + 2 (nw + f sin cos (nx + ga — a) sin ") de 
+ x! cos karx cos + gv — a) sin dx 


n(n7r+q), n det. 


| sin cos (nx + ga — a) cos 


Consequently, since 


sin 
n n 


it follows from (125), (126), and (127) that 

In a similar fashion it may be shown that 

(129) sin kre sin — a) da = ) 


and consequently, since from Theorem I 


(130) d,>Cn, 
*Cf. PICARD, Traité d’ Analyse, 2d ed., vol. 1, p. 256. 
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we have from (124), (128) and (129) 


1 
(131) f sin kv P(X, x) de = YC"). 
Combining (122), (123), (130), and (131) we obtain 

(182) a sin «)dx = 


Since we know that the interval between two successive positive roots of 
equation (107) approaches 7 as a limit, it follows without difficulty that a posi- 
tive constant M can be determined such that 


(185) < Mn. 

Combining (121), (132), and (138) we get 

(134) 


where X, is a positive constant. But since the Bessel function of any order has 
a maximum absolute value for all real positive values of the argument, we have 
from (1052) 

K, 
(135) 


where X, is a positive constant. Combining (134) and (135) we have 


K 


K being a positive constant. Hence the series that we are discussing converge 
absolutely and uniformly in the interval 0 = a = 1. 

Since the developments of the functions (118) in series of the form (119) con- 
verge uniformly in the interval 0 = 2 =1, we have a right to multiply them by 
ax (x) and integrate them term by term from 0 to 1. Performing this opera- 
tion we obtain in view of (117) 


1 
[ sin (2) de = 0 (k=1, 2,3,---). 


Consequently, since any finite and integrable function whose Fourier’s constants 
all vanish is zero at every point at which it is continuous,* it follows that 
x’ (a) and hence x(x) is zero at every point of the interval 0=a2=1 at 
which it is continuous. 


*Cf. HURWITZ, Mathematische Annalen, vol. 57 (1903), p. 440. 


n 
= 
x) | < mi? 
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TueorEM V. Jf f(a) satisfies the conditions of Theorem II and if 
Ay» Ags are the roots, positive or zero, of equation (107) arranged in 
order of increasing magnitude, then the series (106) will be summable to the 
value (110) at every point of the interval 0 <x <1 at which f(x) is continu- 
ous or at which it has a finite jump, provided the latter points lie in an inter- 
val at every other point of which f(a) is continuous; and will be uniformly 
summable to the value f(x) in every closed interval lying in the interval 
0 <2 <1 which does not include a point of discontinuity of f(x). 

We have from equation (111) 


(136) 2)de=0 (k=1, 2, 3, 


where (2) is a function that has the value of the series (106) at every point at 
which /(2) is continuous (see page 421). 

By hypothesis f(a) is finite and integrable in the interval O=a2=1. We 
know from the definition of ¢(a) that (a) is finite and integrable in the interval 
0 <&6=a=1 where 6 may be taken as small as we please. Moreover, it fol- 
lows from Lemma 7 that $( 2) cannot become infinite at the point = 0 to an 
order higher than x~!. Hence $(2) is finite in the interval 0 =a = 1 except 
in the neighborhood of the point « = 0, and a#¢(x) is absolutely integrable 
from 0 to 1. 

Consequently, in view of Theorem IV, it follows from equations (136) that 


(137) x(x) = = 0 

at every point of the interval 0 =«2=1 at which it is continuous. Now we 
know from the definition of ¢(a ) that ¢(2) is continuous at every point of the 
interval 0 <a# <1 at which f(x) is continuous. Hence it follows from equa- 
tion (137) that at every point of the interval 0 <a < 1 at which f/(2) is con- 
tinuous the series (106) is summable to the value (110).* 

We consider next the points at which f(a) has a finite jump and which lie in 
an interval at every other point of whigh /( 2) is continuous. We know that at 
every interior point of that interval, except the point at which /() has a finite 
jump, the series (106) is summable to f(x). Moreover, that series is the sum 
of a series that is uniformly summable and continuous in the neighborhood of 


* Since f(z) —¢(2) at every point at which f(z ) is continuous and since from (109) 


(x) =f(r) + 
n=1 
it follows that 


(2) 

n=1 
is zero at every point at which f(z) is continuous. However, since the series (a) defines a 
function that is continuous throughout the interval 0< z=1, it follows that it defines a func- 
tion that is identically zero in that interval, so that ¢(x) =f(x). 


. 
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that point plus two series that are essentially Fourier developments of functions 
having a finite jump at that point. Hence the finite jump in the function 
defined by the sum of the two Fourier developments has the same magnitude 
as the finite jump in /(2). However, it follows from Frsér’s work that the 
Fourier developments are summable at a point at which the function they define 
has a finite jump, to a value half way between the limiting values of the fune- 
tion. Consequently the series (106) is summable to the value (110) at the 
points which we are discussing. 

That the series is uniformly summable to the value f(a) in every closed 
interval lying in the interval 0 <a <1 which does not include a point of dis- 
continuity of f(a), is an immediate consequence of Theorem II and the fact that 
the series is summable to the value f(x) at every point of such an interval. 


§ 4. Convergence factors and applications. 
TueoreM VI. Jf the series 


(138) u(x) + + + 


is uniformly summable in the interval a = x = b and has the value f(x), then 
the series 

(139) u, (x + + 4 

where Ay, are the successive positive roots of (1), will be absolutely 


convergent in the region 


= 
(140) 4 


and its value u(a, x) will approach f(x) uniformly for all values of x in the 


interval a = x = b, as a approaches + 9. 


We will first show that the convergence factors 
(141) on. ha, 


satisfy the conditions of a general theorem due to Bromwicn.* It will then 
follow that the series (139) is absolutely convergent in the region (140) and that 
its value approaches /( 2) for all values of ~ in the interval @ = x = b as a ap- 
proaches + 9. That this approach is uniform can easily be seen by a slight 
modification of the reasoning in BromwIcn’s theorem. 

BROMWICH’s conditions are 


(a) n|A?v |< K 

a> 0, 
(B) lim nv, = 0 
(y) lim v, = 1, 


a=+0 


~ *Cf. Mathematische Annalen, vol. 65 (1908), p. 351. 
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where w and v are any two integers, v,, v,, v,, ---, are the convergence factors, 
A’v, represents their second differences, and X is a positive constant. That the 
functions (141) satisfy (8) and (y) is at once obvious. It only remains to 
examine (@). 

In order to establish this latter condition we will prove that the second dif- 
ferences of the functions (141) satisfy the inequalities * 


M 
(142) | < Na’ + (O<na=e;n 1,2, 3,---), 
L 
(143) | < (a>0; 3,---). 
I~n a 


It will then follow without difficulty that (a) holds.+ 


Let 
Then by Theorem I 
and hence 
But we have 
(146) Ana + (An +2yn)a a (0<4, <2), 
and since, in view of (144), y, has an upper limit, we obtain 
(147)  DeAntynda 4 e-Ant2yna — 
where JV is a positive constant. 
We have furthermore from the Law of the Mean 
(148) — e—[Ant (gy ( | <1 ), 
from which we obtain with the aid of (145) 
2Kka Ma 
(149) | 2ynla | = 


where V is a positive constant. 


*It may be noted that the inequalities (142) and (143) are analogous to conditions (d) and 
(e) in a theorem on convergence factors proved by the writer (these Transactions, vol. 8 
(1907), p. 300), which isa special case of BRoMWICH’s theorem. The necessity for distinguishing 
between the differences where na =c and na > ¢ often arises in the consideration of convergence 
factors that occur in practice. 

+ The proof of this fact in the present case is analogous to the proof in the case where condi- 
tions (d) and (e) of the theorem mentioned in the previous footnote hold. This latter proof is 
given on pp. 353 and 354 of BRoMwicn’s article referred to above. 
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Finally from (147) and (149) we have 


+ | @—An+ea (Ant | Na? Mx 
+ 


and hence the inequality (142) is satisfied. 
We obtain from (146), (148), and (145) 


—An —An+ = Ant 2 —[Ant 2ynt \ 


a? 2Ka _ 2y? 2K 
(A, + + A, + 27, + 8,( Wa + a 
2K 
L 


where Z is a positive constant, and therefore the inequality (143) holds. 

Consequently, as we have pointed out above, our theorem is proved. 

An analogous theorem about slightly different convergence factors is the 
following : 

TuHeoreM VII. Jf the series (138) is uniformly summable in the interval 
a =2=b and has the value f(x), then the series 


242 \ 223 
(150) 4+ w,(a) + 4 
where X,, r,, %,, +++ are the successive positive roots of (1), will be absolutely 
convergent in the region (140) and its value u(a, x) will approach f(a) uni- 
formly for all values of x in the interval a = x = b, as a approaches + 9. 
This theorem also can be proved by showing that the convergence factors 


satisfy the conditions of Bromwicn’s theorem. The proof is quite analogous to 
the previous proof, save that instead of the inequalities (142) and (148) we 
establish the slightly different inequalities 


(152) | A? | < Na? a—1, 9, 3, 
(153) |< (a>0; 


The results just obtained have a direct application to many problems in mathe- 
matical physics. 

Problem I. The problem of determining the flow of heat in an infinite, cir- 
cular cylinder, when the initial and boundary conditions depend only on the 


‘ 
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perpendicular distance from the axis of the cylinder, reduces to the determina- 
tion of a continuous function * of r and z, v(r, z) such that 


154 Ov Cv 1 Ov 
(154) \ OF + 
in the region 

5B z>0 
(155) 
and 

Ov 

(156) t+hv=0 (z>0,r=1), 
(157) v= f(r) (z=0), 


where a, /, and / are constants and f(r) isa function of r that gives the initial 
distribution of the temperature. 
It can be shown easily that the function 


(158) F(X, 

where A is a constant and X is a root of the equation 

satisfies conditions (154) and (156). It remains to construct from functions of 
the form (158) a continuous function which will also satisfy (154) and (156) 
and in addition will approach continuously the boundary value /(7), where /(7) 


is continuous, and which will remain finite where that is not the case. 
Consider the series 

(160) 7) + 7) + 
where the d’s are the roots, positive or zero, of equation (107) for v= 0 
arranged in order of increasing magnitude and the A’s are the coefficients of 
the series of the form (106) corresponding to the given function f(r). If 
the function f(r) is finite and integrable in the interval 0 = 7 =1, we obtain 
from (46), (50), and (105a) . 

|A, Fi(A,,7)| (0=r=}), 
where (A is a positive constant. Hence for the general term of (160) we have 
in view of the first inequality in (23) 


* By a continuous function we mean that v(7r, z) is to be continuous in the two variables r 
and z throughout the region 
z>0 
(4 =r= 1) 


and is to approach the boundary value f(r ) continuously wherever f(r) is continuous. 
t Cf. CARSLAW, Introduction to the Theory of Fourier’s Series and Integrals, etc., Chapter XV. 
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Kv x, K Kk K ( ) 


| A e—aazz F ( r r) | 


an? arniz ~aCiniz” aCiz,ni 


and therefore the series (160) is uniformly convergent in the region 


(161) 


Since, moreover, the terms of (160) are continuous functions of + and z for all 
values of r and z, it follows that the series defines a continuous function of r and 
z throughout the region (161), and therefore, since z, is any constant greater 
than zero, throughout the region 


0=r=1 
(162) ( met ). 


In similar fashion it may be shown that the function defined by (160) satis- 
fies conditions (154) and (156). It only remains to consider the behavior of 
the function as we approach the line z = 0 from the positive side. 

We know from Theorem V that if (7) is a function that is finite and integ- 
rable in the interval 0 = x = 1, then the series 


(163) A, 7) + A,Fy(Ags 7) 


where the A’s and )’s are those of the series (160), will be uniformly summabie 
to the value of f(7) throughout every interval that lies within a larger interval 
in which f(7) is continuous. It follows then from Theorem VII} that, as z 
approaches zero through positive values, the function defined by (160) will 
approach the value f(7) uniformly throughout any r-interval lying in a larger 
interval in which f(7) is continuous. Consequently, if 7’ is any interior point 
of an interval in which /(7) is continuous, the function defined by (160) 
approaches the value /(7) as we approach the point (r = r’, z = 0) along any 
path lying in the region (155). 

We now wish to examine the behavior of the function defined by (160) when 
we approach a point that is not an interior point of an interval of continuity of 
J(7). If we define S\(7) for the series (163) in the same way that S (2) 
was defined in (34) for a corresponding series, we may, for values of z > 0, 
write the series (160) in the form ¢ 


n=1 


*Tf 7,0 we must take the (n+ 1)st term as the general term instead of the nth term in 
order that this inequality may hold. 

{ In this theorem we did not consider the possibility of 4, being equal to zero. However, 
this is an unessential point, as the behavior of a finite number of convergence factors at the 
beginning of the series does not matter, provided they approach 1 as the parameter of which they 
are functions approaches zero. 

¢ Cf. equation (7) of the writer’s article in these Transactions referred to on page 430. 
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If now we make the change of variable 


the expression (164) reduces to 


(165) 2S, (7) — 4 ], 


n=1 
But we have from the corollary of Theorem II * 
|S(r)| <M (n=1, 2, 3, ---; 
where M is a positive constant. Consequently, since the convergence factors 


(151) satisfy the conditions of Bromwicu’s theorem, the series (165) is, through- 
out the region 


less in absolute value than 


n | — 4 | < MK, 


n=1 
where A is a positive constant. Hence if we approach any point (r= 7’, t= 0) 
such that 
1, 
along any path lying in the region (155), the function defined by (160) remains 


finite. 
Problem II. The problem of determining the steady flow of heat in a semi- 


infinite circular cylinder whose base is kept at a temperature 

v= f(r), 
where r is the perpendicular distance from the axis of the cylinder, reduces to 
the determination of a continuous function + of r and z which satisfies 


1 oO Ov 


in the region (155), and also satisfies (156) and (157). 


It can be shown easily that the function 
Ac™ F(A, 1); 


where A is a constant and 2 is a root of equation (159), satisfies (166) and 
(156). We form now the series 


* This corollary was proved with regard to the series (35), but as (35) differs from (106) by 
at most one term that is finite for all values of the variable involved, it holds for the latter 
series also. 

+ Cf. note, p. 432. 

¢ Cf. CARSLAW, loc. cit.. p. 321 


a>0d 
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0 


where the X’s are the roots, positive or zero, of equation (107) for »=0 
arranged in increasing order of magnitude, and the A’s are the coefficients of 
the series of the form (106) corresponding to the given function /(7). 

If f(*) is finite and integrable in the interval 0 = r = 1, the series (167) will 
converge for z > 0 to a continuous function of 7 and z that satisfies (166) and 
(156), and approaches the proper boundary value when we approach a point on 
the line z = 0 that lies within an interval of continuity of f(7). Furthermore 
the series (167) remains finite when we approach any point on the line z = 0 that 
lies in the interval 

0<cH=ri=l. 

The proof of the above facts is analogous to the proof in the previous prot- 

lem, except that we use Theorem VI instead of Theorem VII. 


THE UNIVERSITY OF CINCINNATI, 
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SINGULAR POINTS OF ORDINARY LINEAR 
DIFFERENTIAL EQUATIONS* 


BY 


GEORGE D. BIRKHOFF 
Let 


(1) dy, 


j=1,2, 


be a system of n ordinary linear differential equations of the first order in which 
we have for |x| > R 


1 
(2) (x) = a,x" + + + af (8, f7=1,2,--+, n), 


q being an integer. We propose to investigate the nature of the solutions of 
(1) in the vicinity of a = 0, which is taken to be a singular point of the given 
system of equations so thatg =—1. The integer g + 1 is termed the rank of 
(1) at # = oo, after Porncaré.f It will be further assumed that the n roots 


As a,, 
of the characteristic equation 
(3) — = 0 
[where 6,, = 0 (i +), and 6,,=1] are distinct, and in the case g = —1 do not 


differ by integers. The writer hopes to discuss the general case later. 
It is apparent that any proper linear transformation 


(4) 


in which the functions A,,(a) are analytic or have a pole at x = 0, will trans- 
form (1) into a system of equations (1) with coefficients @,,(«) of the same form 
(2), although ¢ is not necessarily equal to g. Any such system of equations (1) 
will be said to be equivalent to (1) at x=. 

The first part of the paper deals with the determination of the simplest sys- 
tem of differential equations which is equivalent to (1) at ~=oc. In $1 an 
important Lemma on Analytic Functions is stated and proved. The applica- 
* Presented to the Society (Chicago), April 17, 1908, under a different title. 

tActa Mathematica, vol. 8 (1886), p. 305. 
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tion of this lemma in § 2 shows that there always exists a canonical system of 
equations equivalent to (1) at 2 = o having the form 


dy, | 
= Py (¢=1, 2, --+, 2), 


in which P(x) are polynomials. In §3 it is shown that the rank of this 
canonical system of equations at 2 = oo can be taken to be qg so that the degree 
of the polynomials P,,(x) does not exceed g + 1. 

The significance of these results is best seen by considering the case of a reg- 


ular singular point g = —1. The canonical system of equations equivalent to 
(1) is then in simplest form (§ 3) 
dY, d¥, _aY, 


and a fundamental system of solutions is 


7(2 2 
Y® Y® == Vi 


. . . . . . . 


Now the solutions of (1) and the above system of equations are connected by a 
relation 


DA, (2)¥, (t=1,2,---,m), 
j=! 


in which the A,,(a) are analytic at 2 =o (or have a pole, but this case is 
easily excluded). It follows that there exists a fundamental system of solutions 
of (1), (= 1, 2, ---, m), of the form 


( LY =r, (x) (é,j—=1, %). 
k=1 


This is the fundamental theorem for a regular singular point. The solutions 
of the canonical system of equations equivalent to (1) at x = c play the same 
réle in the case of an irregular singular point (q = 0) that they do in the case 
above considered. 

The nature of the solutions of the canonical system of equations can be dis- 
cussed by special and elegant methods when g=0. In the case g=0 it is 
known that the solutions of any system (1) with rational coefficients a,,(a) can 
be expressed by a Laplace integral 


J en,(z)dz (i=1,2,---,m), 


n 
n 
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and that from this formula can be deduced the behavior of the solutions in the 
vicinity of « = oo. This method of treatment is due to PoIncaré,* and has 
been developed more fully by Horn.} It is applicable in an extremely simple 
form to canonical systems of equations for which g=0. If g > 0, generalized 
Laplace integrals of a simple character will be shown to be available. The 
integrals for g = 0 are considered in § 4. 

These integrals are employed in §5 to find the asymptotic nature of the solu- 
tions of the canonical system of equations equivalent to (1) at 2 = oo in certain 
sectors of the a-plane. The asymptotic form of the solutions of (1) in the same 
sectors follows as an immediate consequence. The powerful and general method 
of successive approximation which Horn { has employed in various cases will no 
doubt lead to the same results. 

A Riemannian characterization of the solutions of the canonical system of 
equations in the vicinity of « = oo is given and extended in § 6 to the solutions 
of (1). The number of characteristic constants in the solutions is found to be 
the same as the number of arbitrary constants of the canonical system of 
equations. 

The characterization suggests the probable generalization of the Riemannian 
problem when the singular points of the system of differential equations are not 
regular. This formulation and a count of constants is given in § 7. 


§1. Auziliary lemma. 


In order to obtain the simplest possible system (1) equivalent to (1), we can 
choose as we please the n’ functions d,, of the transformation (4) but with the 
restriction that they are analytic or have poles at 2 = oo and that the determi- 
nant |A,,| is not identically zero. How can one choose the functions 2,, so as to 
free as far as possible the solutions of (1) from singularities in the finite plane ? 
The answer to this question, which will be fundamental in our study of (1), is 
derived from the following lemma. 

LEMMA ON ANALYTIC FUNCTIONS. Let 0,,(x) be a set of n° functions of x, 
single-valued and analytic for |x| > R but not necessarily analytic at x = co, 
and also such that the determinant 


+ 0 (|z|>2). 


Then it is possible to so choose n* multipliers (x), analytic at x = co, that 


> (2) = (x), (or =="E,(2)) 
(i, k=1,2,---,m), 


*American Journal of Mathematics, vol. 7 (1885), pp. 203-258. A single linear 
differential equation of the nth order is treated in this and other papers referred to later. 

t Mathematische Annalen, vol. 50 (1897), pp. 525-556. 

tCrelle’s Journal, vol. 133 (1907), pp. 19-67. 
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where m is a fixed positive integer or zero, and the —,.(x) are entire functions 
of x for which the determinant |&,(0)| + 0. 

Proof. The two alternative forms of the Lemma are clearly equivalent. 
We shall confine ourselves to the proof of the second of these, which is the 
easier to give in matrix notation. 

In order to elucidate the lemma let us consider first the case n=1. We 
have to determine i,, so that 


Since @,,(x) + 0 for |x| > the function log is expansible in a Laurent 


series for |x| > R except for a term a log x, 


Te 
log = >> + alog x. 


v=-D 
Take then 
- > > 
ine m= a =e"? 
if a<0,or 
+2 
A, =e m=0, = ate 


if a = 0, and we have the desired solution. This explicit treatment is not pos- 
sible if n> 1, in which case we shall have to construct a solution by an infinite 
process. 

As it is convenient to use matrix notation, let @ denote the matrix || 6, || and 
@~' its inverse matrix ||0>'|| which is of the same form, since |@,,| + 0 for 
|2|> Also we shall write A = 

Furthermore, if ¢ is single-valued and analytic in the vicinity of « = oo , so 
that it is expansible in a Laurent series, 


= a,x", 
we write 
+n 
P($) = 24,2", N(¢)= 4,2", N,,(¢)= a, 2”. 


Here m is some positive integer whose value will be specified later. We have 


the obvious relations 
$= P($)+ W($)=P($) + (6), 
P(N($)] =P[N,(¢)] = NL P($)] = 
More generally, if we have a matrix = ||¢,||, we shall use the notation 


Trans. Am. Math. Soc. 29 
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Let now 7'=||¢,,|| be any matrix of constants of determinant not zero. 


Consider the matrix equations 
(5) r=P(e"'-V) 


exists, we will have at 


If a solution of these equations T = ||¥,, 
once the desired matrix 


(6) A=N(@"- WV). 


To verify this statement we notice first that the functions X,, are analytic at 


x = 00 by (6), since they are convergent series of negative powers of x. More- 
over from the first equation of (5) we obtain 


(7) r+ 
From this equation (7) we derive 
(8) 


on multiplying by the matrix @. Furthermore, by the second equation of (5), 
we obtain also 


(9) 
in which the elements of # are polynomial of degree m —1 at most in 1/2. 
If we add (8) and (9), we find by means of (6) 

6-A = — I) — Rf. 


The right hand member forms a matrix ||2~"€,(a)|| in which the &,(a) are 
entire functions of « with the desired property 


|E,(9)| = |¢,| +9. 


We have thus proved our statement. 
It remains to construct a solution of equations (5) for some m. This we 
do by a method of successive approximation. 


Write 
= Tz, ° 
r,= P(e"-¥,), 
(10) V,= + N_(98-T,), 


= 
Tx" + N_(9-T,), 


If the sequence of matrices Vv: r tend to limit matrices, it will not be difficult 
to prove that these limits form the desired solutions. 
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Define the difference matrices 


r,=",, 


These matrices for p > 0 satisfy the equations 
(11) T.=P(e"-¥). 


It follows from the first of these equations that VY contains only negative powers 
of x. Hence the positive powers which come from the multiplication of @~' and 
Vin the second of these equations arise from the positive powers in the ele- 


ments of @~' and from the negative in VY) so that 


(12) P(@"-¥ )=P[P(@")-¥ 
Likewise we find 
(13) V,,(@-T,_,)= 


since the elements of I, , contain only positive powers of x and a constant term. 
Thus (11) can be written 


(14) v= r= P(@")-V] 
In this form we shall prove the convergence of 


for large enough m. 
Choose r,, /2, so that 
then we can choose JM so large that* 
M 
1_” 


0 
since P(6;') is an entire function of x. Moreover since V(@,,) is analytic for 
> &, the series for V(@,,), 
a) 


| a 


* The notation, a< 3, means that the coefficients of the various powers of xz in the power 
series a do not exceed the corresponding terms of 3 in absolute value. I shall extend this nota- 
tion of POINCARE to matrices A and B of which corresponding terms aij and {4 are power series. 
Thus the notation A < B implies that aij < $i for all values of i and j. 
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is convergent for |~| > 2, so that 


lim =0. 


p=2 0 


Hence for large enough m we have | p> m, no matter how small 


e be taken. Thus we can choose m so large that for any small e 
Take m so large that 


(17) 


It is to be noticed that both the comparison series used converge for 
r,>|z|> 
Now choose Z so large that 


(18) r=T,« 


First by (16) and (18) we see that 


neLr, R, , neLr, R, 1 1 


| neLr, R 
Ry) (x — hy) 
Next, by (15) and the above inequality, we see that 


R, | n'eMLr* R, 


as before. And thus in general 


as may be proved by induction. If (17) holds, we shall have convergent series 
for y, and ¥,, if 
Furthermore, from the equations (10) we infer that 


r= P(e"), 


* The compound matrix has n terms in each element ; hence we get a multiplier n. 


L 
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so that y,, must be an entire function of x. Then likewise we infer that 
Thus || 7,,|| and || ¥,,|| satisfy (5). 


§ 2. Hwxistence of an equivalent canonical system. 


We are now in a position to answer the question which was raised at the 
beginning of § 1. 

By the hypothesis concerning (2) there are no finite singularities of the func- 
tions a,,(x) of (1) exterior to the circle |x| = #. Let x describe a closed path 


exterior to this circle about # =o. It is easy to show that there will exist n 
linearly independent particular solutions 


of’, 
which undergo the substitution 


(n) 


unless two or more values of p are equal, when degenerate cases may occur.* 
We will suppose for the present that there exists a fundamental system of 
this character (19) for which, of course, 


Define now functions 0,,(a) by the n’ equations 
= 28, (2), (4,= log ) 
2rV—1 
It may be inferred at once that the 0, are single valued and analytic for 


> # and that 


| 6 (x) | emit +hy) oS [4 +A 


Thus the 0,(2) fulfill the restrictions of the Lemma. Hence we can deter- 
mine a set of functions ,,(«) analytic at 2 = oo in such a manner that 


the being entire functions of and |&,(0)| + 0. 


* Tt is here simply a matter of the canonical form of a linear substitution. 
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The transformation (4) changes (1) into a new system of equations (1) with a 
fundamental system of solutions 


Therefore we can transform a system (1) into an equivalent system whose 
solutions have no other singularity in the finite x-plane except x = 0, and this 
singularity will be of the above simple form. 

The new system of equations can be written 


dj, _ 
de 
(20) dx In 0 (i=1,2,---,n). 
i | 


Hence we see that the coefficients d;,(«) are quotients of entire functions, and 
therefore meromorphic in the finite plane with poles only at points for which 


A= = 0. 
On the other hand, by a fundamental property of any transformation (4), the 
new system of equations (20) is of the same form as (1) so that the @,(a) are 
analytic at # = oo, or have a pole at that point.* 
These two properties of the functions @,, show that they are rational functions 


of x. Thus the new system of equations equivalent to (1) at x= 0 can be 
written 


dy, n 
(21) LR, 
in which the coefficients R(x) are rational. Furthermore it possesses a funda- 
mental system of solutions ° 


(22) = E.,( ) 
where the £,(x) are entire functions of x and |£(0)| + 0, and the poles of 
f(x) lie at points for which 

A=|y?|=0. 


* Write the transformation inverse to (4) in the form 


in which, of course, the 4;(xz) are analytic at x= 0, or have a pole at that point. A direct 
computation shows that we have 


n 
ay= 
_= 


Ain dn Ay + Ang Aik. 
dx 


k 1 


Yis 
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By a succession of transformations (4) we shall now remove all the singular 
points of (21), not at e=0 and x=oo. These can only be apparently 
singular points because of the character of the fundamental system of solu- 
tions (22). 

Since from (21) we have 


ced + Ryn (x) 


it is clear that the entire function A is the product of a factor which never van- 
ishes, and of factors of the form (#—~2’)', X being a positive integer. The 
singular points will lie at x =~’. These last factors can be removed as follows. 
Choose constants c,, ¢,, ---, ¢,, not all zero, so that for x = a’ 


(23) +---+e,y=0 (j=1,2,---,n) 


This is possible because by hypothesis A= 0 forx=a2’. Let c, be one of the 
quantities c; which is not zero. Make the further transformation of the form (4) 


Yi= Yi (i+k), 


The fundamental system of solutions 7/) of (21) is transformed into another 
such system of the same character, 


(i+#k), 
j= + + - 6,9) 


; 


Accordingly the new system of equations has again the form (21). But we 


have 
im 
so that a factor x — x’ has been removed from A. 
In this way one can remove in succession each factor x — 2’, and finally reach 
the case when A does not vanish except perhaps for2=0. Thus a system of 
equations (21) exists whose solutions have the additional property 


(24) + 0 (240). 


It follows that the coefficients /?,,(#) cannot have any poles except at x = 0 
and 2 = oo, and so are either polynomials or polynomials multiplied by a posi- 
tive integral power of 1/2. The index of this power does not exceed 1. For 
suppose we have 


1 T 
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in which 7 > 1 and not every ¢, is zero. Substituting ¥/) = x“, (a) in (21) 
and equating coefficients of lowest powers of x, we obtain 


0= 2 (i, 2,---,n). 


By taking i = 1, 2, .--, m successively it is inferred that ¢, = 0 (i, 7 = 1, 2, - --, n) 
since |&.(0)|+ 0. This is contrary to hypothesis. Therefore we must have 
y 


P;;(x) 


where P(x) denotes some polynomial in . 
THeoreM. There always exists a canonical system of equations equivalent 
to (1) atxr=o, 


dy. 
(25) = > P,(x)y, (i=1, 3, +--+, 0), 


in which the P,,(x) are polynomials. 

The proof of this theorem has been based on the assumption that there exist 
n linearly independent solutions of (1) of the character (19). The extension 
of the above proof to the degenerate cases is illustrated by taking a degenerate 
ease for which p,_,=p,. Then there exist n linearly independent solutions 


such that we have in place of (19) 


Yo = PiY2's = py, 


n—1) 


n—1) —~(n—1) n—1) (n—1) 

Yi Y2 Yu = Pn-1Yn 
n—1) — n) n—1) n n—1) 
In this event we may write 


= 28, (x), (i=1, 2,---, j=1, 2, n—1) 


= [ (2) + 6, n—1(#) log | (i=1, 


—1 Pn—1 
and the functions @,(2) thus defined will be single-valued and analytic for 
|e|> R. 


Furthermore we have 


| 


Thus the functions ,,(a) satisfy the restrictions of the lemma of § 1, and a set 
of multipliers A, (x) will be determined. If now the transformation (4) be 


n 
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made, the new system of equations in 7, will have a fundamental system of 


solutions 
oP = E(x) (i 1,3, 


1 
= E(x)+ E, (i <->, ah. 
2rp, 


n 


in which £,,(2) is an entire function and | &(0)| + 0. The rest of the argu- 
ment in the degenerate case is very like that employed before. 


§ 3. Formal solutions and the polynomials P,(x). 


It can be shown easily by means of the formal solutions of (1)* that the poly- 
nomials P;,(«) cannot all be of degree less than g + 1. 

Let us first consider these formal solutions in the case g = 0. If we substi- 
tute in (1) formal expressions 


y, == A, (2) (i 1,2,--*,%)s 
itl 


1 
A ,(2) = A, + A+ 


and attempt to determine a, 8, ---,r,«, A,, A’, --- so that the equations 
resulting are formally satisfied, and not all the » quantities A, are zero, we 
shall find that a is one of the roots of the characteristic equation (3) and that 
we can then determine 

A, B, AM, 
in succession ; thus there arise precisely n formal solutions of (1) 


where 
SP = eu A, (x) (i,j =1, 2,---,#), 
mitt 


1 
A,(x)=A,+ A} 


ion 


The first set of equations determining A,, are 


(27) A,a,= (i, j=1, 2,---,m). 
k=1 


* SCHLESINGER, Vorlesungen iiber lineare Differentialgleichungen, pp. 176-182. For a discus- 
sion of regular singular points see pp. 122-171. 
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It is necessary for us to know also that the formal determinant 


1 


does not vanish identically. As a matter of fact, the constant D = | A,,| is not 
zero. Otherwise we could find multipliers m,, m,,---, m,, not all zero, such 
that 

n 

= (i=1, 2,---,m). 


Multiplying equations (27) by m,(j7 = 1, 2, ---, m) and adding, we obtain 


n 


A,,4,m, = 0 (s=1, 2,---,). 
j=l 
A second set of multipliers is = likewise a third set is 


Thus we obtain » sets of multipliers My a" It is clear the 
that we have 


A,, m,[¢, =0 

for all values of ¢,, ¢,,-+-,¢ 
that 


Since a,, +++, are distinct we conclude 


A,m,=0. 
Hence for any value of j for which m, + 0 we have 


A,,=0 (i=1, 2,---,n), 
which is impossible. 
Let us turn now to consider the formal solutions of any system (1), equivalent 
to (1) at 2 = co by means of a transformation (4). The formal solutions of (1) 
are defined by the equations 


SP SY (i, j=1, 2, ---,n). 


Hence S‘/ contains the same exponential factor as SY), although x, may differ 
from «, by an integer. The formal determinant | SY?| is also not identically 
zero, since we have 
[SP] = 877 


The same results are obtained for the case g =— 1 except that the exponential 
multipliers e” are lacking in the formal expressions SY) and S‘). 


We will now prove for the case g = 0 that the rank of any system of equa- 
tions (1) equivalent to (1) at «= 00 is at least g+1. Otherwise we must 


n 
n 
k=1 
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have equations like (27) (a‘') being replaced by @;')) to determine the constants 
A,, of the formal solutions of (1), in which however 


a..=0 f=1, 3, --*, 
These equations are 
A,a,=0 (i, j=1,2,---,m), 


whence we see that A, = 0 (i=1,2,---,n) for any a not zero. This is 
impossible. 

In particular, then, the rank of the canonical system (25) equivalent to (1) at 
x= oo isat leastg+1. Zhe polynomials cannot all be of degree less 
thang +1. 

We shall now prove that the degree of the polynomials P,,(x) can be taken 
not to exceed q+ 1. 

Let us denote the » formal solutions of (25) by ZY), TY’, ---, TY 
(j=1,2,---, where 


TY = (a), 


(28) 
= B,, BY) - 
42 


x,—J, being an integer. Let ¢; be the index'of the highest power of 1/x 
which all the B(x) (j =1, 2, ---, m) contain as a factor for fixed i. Then 
we may write 


1 


io 


By a succession of transformations (4) of such a character as not to change 
the form of (25), we now proceed, if necessary, so to transform (25) that 


(29) +0. 
Assume the indices 1 to n so chosen that 
=t,=---St,. 
If |C,,| = 9, one can choose constants c,, ¢,, ---, ¢, not all zero, such that 
¢,C,, +¢,C,, + =0 (j=1, n). 


Let c, be the last one of this set which is not zero. The transformation of the 
form (4), 
y (i+s), 


a t,—te 
Y,= 'Y, + +CY,> 
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leaves (25) unaltered in form. The new formal expressions 7'/) are the same 
as 7’) fori+s. Also 


8 1 

i=0 


The terms (7 =1, 2,---,) are all zero by hypothesis. The for- 
mal expressions TS (j =1, 2, ---, m) will contain another factor 1/x, and we 
have increased ¢, by unity or more. 


But we have 
| =<, 


directly from the transformation. In this way we can increase ¢,, t,,---, ¢, 


without altering, save by a constant factor, the formal determinant |7’))|. 
Since the determinant contains 1/2 at least to the power 


we see that this process of transformation must finally come to an end. This 
occurs when (29) is satisfied. 

When (25) has been transformed so that also (29) holds we may prove at 
once that the degree d,, of P,, can not exceed g+1+t,— t,. In fact, if some 
of the integers 

d,—(q+1+t,—¢) (i, j=1, 2,---, ») 
are positive, let t be the greatest of them. Also denote the coefficients of the 
term 
in P(x) by p,. Substitute the formal solutions in (25). When coefficients 
of highest powers of x are equated to zero one obtains 


Cp (4, 9, 8) 


which is not possible in view of (29). 

We note at this point that the degree of P;,(x) does not exceed ¢ +1 if 
i=j. We will now make further transformations so that this will also be true 
when i 

Let us first attempt to remove the constant terms of P,,(x) (j=1, 2,---,n—1) 
by a transformation 


Y; (i n), Yn = Can Yn + Cy 


n 
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where ¢,,,=1. We will then obtain 


P,, (x) = Py(x) —¢,,P,,(«) (i,j#n), 

P,,(2) = x ) (i+in), 

j=l 


P(x) = Pin 
j= 


It is to be noted that the degree of 7. (x) (i,j =1, 2, ---, nm) does not exceed 
qt+1+t,—t. If p?) denotes the constant term of P,; (x) we have to satisfy 
the conditions 


0) (0) — 
Pon Pin = 0 (kK=1, 2, n—1 ). 
I= 


Or if we write 


these conditions may be stated in the form 


Pye — Muy = O (kK=1,2,---m). 
j= 
This gives the equations for ), 
| — 5, =Q. 
ij 


Now each transformation employed since we obtained (21) has not affected 
the form (22) of the solutions. Hence (25) still possesses a fundamental system 
of solutions of this form. A direct substitution of these solutions in (25) will 
show that the roots » are the quantities 


k,—m, +--+, k,—m. 


The roots are therefore distinct and do not differ by integers in the case 
in which (1) has x solutions of the form (19) where p,, p,, ---, p, are distinct, 
inasmuch as 
1 
= log p, (j=1,2,---,m). 

We confine ourselves to this simplest case. 

Choose for \ any root &, — m, determining a set of values Cy (j=1,2,---,n) 
in which ¢ +0. Such a root exists, since otherwise the linear equations 


n—l 
— Ae, = 0 (k=1, 2,---,n—1) 
j=! 


jk 


could be satisfied for m values of X, and there are only n —1 such values. 
The corresponding transformation (c,,= 1) will remove the desired terms. 
=k, — m in the new system of equations. 


(0) 
un 


Also we have p 


nm 
n n 
n 
> (0 
njf jn nn 
n 
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Next we remove the constant terms of P,_, (x) (j7=1, 2,---,n—2) by 
a transformation 


This can always be accomplished as before and we shall have reduced the con- 
stant term of p,_,,_,(“) tok,_,—m. Continuing thus, we finally remove all 
the constant terms of P(x) (i <j) leaving 


(0 « (0) 
PL=k,—m, pR=k,—m, ---, p=k,—m. 


When the constant terms are removed, we apply a transformation 


((+2), Ya + CX" Ys 


to remove any term of P,,(2) (a> 8) of degree less than ¢,—¢,. A 
direct computation shows that we have 


= P,(z) (i+a,j+3), 

= Pig — ca" P,, (2) (i+a), 
(31) 

= — ca" Pg; (a) (#3), 


Pp, ( Fa (x) — [ (x) + cu" (x )]. 


It is clear that the degree of | A (x) (i,j =1, 2,---, m) is not greater than 
q+1+t,—t, since tr ~< t, —t,. 

Let us take 7=1. All the new terms in P, (x) (i<j) are at least of first 
degree, since a multiplier x appears throughout in them. Furthermore, since 
P,,(@)(i>a)and P,;(x) (S8> 7) contain no constant terms after the previous 
transformation, it is clear that there are no new terms of first degree except 
possibly in 


(ix<a), (8<j), Pig (2). 


Also the terms of the first degree in P,,(x) can be made zero by a proper 
choice of c, namely such that 


Pils —¢(k, — kg) =9. 


This shows that we may remove the terms of first degree in Py (x) (t,— ‘> 1), 
by taking them in proper order, namely with increasing i and decreasing j. 

In the same way we may continue, taking t= 2, 3,---, until all terms of 
P.,(x) (i>j) of degree less than ¢;—¢, are removed. We may therefore 


write finally 


where Q(x) isa polynomial fori >. This holds likewise for i =j, since then 
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t,—t,=0. Also Q(x) will be of degree g +1 or less, since P,(x) is of 
degree gq +1+t,— t, or less. 
The transformation 


Y, = (i =], 2,-- n) 
will give us a system of equations 
dy, 
= P..(x)y,, 


in which 
P(x) = P, =Q,,(2) (i+7), P(x) = P,(x) +t, =Q, (x) + t,. 


Thus there exists a canonical system (25), equivalent to (1) at 2 =00, in which 
the degree of P;,(x) does not exceed g + 1. 
We still have at our disposal a linear transformation with constant coefficients, 


F,= y;- 
In view of the nature of the roots of the characteristic equation (3) we may 


simplify P,,(«) by this means so that the equivalent canonical system of equa- 
tions becomes * 


P,(x)¥,, 
dx j=l J 
(32) Pi = + + (i+), 


P(x) p® + Pra + pat + 


The result which we have obtained is that there exists a canonical system of 
equations of the form (32) equivalent to (1) at x=. This form will be 


hereafter assumed. 


In the case g = — 1 the canonical system reduces to 
d¥, d¥, dY, 


But this system of equations is equivalent to (1). Hence, as stated in the intro- 
duction, there exists a fundamental system of solutions of (1) 


where A, (a) are analytic at 2 = oo or have a pole at that point. As there are 
precisely n formal solutions, 


y, = A, (x), Y, = 2% A, (x), = A,,(x), (fj =1, 2,--+, 2), 


1 
A,,(2) =A, + 


* For an analogous reduction in the case g — 0 see SCHLESINGER, loc. cit., pp. 183, 184. 


| 


454 G. D. BIRKHOFF: SINGULAR POINTS OF [October 


we infer that A,,(x) = A, (a) and is therefore analytic at x=0o. This is 
the fundamental result for a regular singular point and is due to SAUVAGE.* 
§4. Generalized Laplace integrals. 


The problem of the paper now resolves itself into a study of the solutions of 
(32) forg = 0. 
Let us consider first the case g =. The equations (32) are then 


= (p® + a,2)¥, + +---+ pF, 


= + pe +--+ 


The formal solutions will have the form (7, TY, ---, 7%) (j=1,2,---,n), 


where 


T = B(x), 
1 
B,(x) = B,+ Bi) 


Direct substitution shows that 
1, = (j=1,2,-++,n) 


In this system of linear differential equations write 
(84) = fen(z)dz (i=1,2,--+,m). 
By employing the formula 
ferx u(z)dz = [e*u(z) ] — fe*du (z), 


in which the path of integration will be subsequently so chosen that the 
bracketed term may be omitted, we obtain, by equating the coefficients of e* in 
the integrands on both sides of (38), 


dn 
dn ( (0) 


dn 


* Annales del’ Ecole Normale, series 3, vol. 3 (1886), p. 392. The analogous result 
for a single equation of the nth order is due to Fucus. 


. . . . . . . . . . . . 
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We shall show how to employ (35), thus formally deduced, to obtain an actual 
representation (34) of the solutions. 

The system (35) has + 1 regular singular points at z=, ---, 
and z = o, and no other singular points. The exponents at « = 2, are all zero 
except one which has the value — p) -1= —/,—1. These results show that 
we may choose a fundamental system of n solutions, of which x — 1 are regular 


at 2 = a, and one has the form 
7% = (2 (2 = a; ) 


in which [/)(z — a,) is regular atz=a,. This is true except possibly when 
J, is an integer. We will refer to this possibility later. 
The formula (34) will represent a solution of the given canonical system of 


equations provided we write 


and the corresponding contour c, of integration is taken to be a “loop circuit” 
that is, a cireuit extending along a (properly chosen) ray from z = oo to the 
vicinity of z= ,, and then around the point in a negative sense, returning 
finally to z = co along the ray. The direction of the ray must be so chosen 
that it does not contain another singular point z = a,, and that 2?(«z) (the 
real part of xz) is negative along the ray. When this is done there is no diffi- 
culty in proving that each set of formulas 


represents a solution of (33), at least in certain sectors of the plane. 

In the case in which — 7, — 1 is a positive integer the loop circuit is replaced 
by a straight line path from z= 0 toz= 00. When —/;,—1 is negative or 
zero, there will be a solution which contains a logarithm. We choose this solu- 
tion to replace 7, nf, ---, 7. In all cases we obtain a formula for each 
singular point. 

This elegant method of attack (for a single linear differential equation of the 
nth order with rational coefficients) was given by PorIncar&é in a fundamental 
paper.* The method has been more fully developed by Horn,} who shows that 
a direct study of the solutions can be made by these formulas in the complete 
vicinity of 2 = oo. 

When g > 0 the above method fails. Porncaré has indicated a method of 


* American Journal of Mathematics, vol. 7 (1885), pp. 217-232. 
{ Mathematische Annalen, vol. 50 (1897), pp. 525-556. 


Trans. Am. Math. Soc. 30 
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quadratures reducing the case g > 0 to the case g=90, but the application 
of this method seems impracticable. * 

It is the aim of this paragraph to show that, at least in the case before us, 
formulas exist for g > 0, analogous to those given when g = 0.+ 

First let us consider the case g = 1. The equations (32) are in this case 


= + pile + 20°) + (Py + pe) phe), 
dY, (0) 1) (0) (1) 2\ (0) (1) 
(36) dx ~ (Pri + Pi”) 1+ + + (Pp 2n + pi”) 
d¥,, 


d: 


The formal solutions are ( ---, (j =1, 2, ---, n), where 


x2 


B,,(x) B,, + B') + 


If these formal expressions be substituted in (86), we obtain 
(37) B, = p™ (Sunt, 0). 


We shall take as the generalized form of (34) in this case 


(38) [ng(2) + (2) (i=1, 2,--+,n). 


Now substitute this expression (38) in (36), integrating by parts according to 


the formula 


in order to reduce the terms which contain x* and 2’ to terms in z' and 2° re- 
spectively. If the bracketed terms are omitted and the coefficients of e”** and xe** 
in the integrands are then equated, there result the equations : 


dn, 
— (22 —a,) na + PY) (¢=1, 2, ---, n), 


dz 


* HorN, Acta Mathematica, vol. 23 (1900), pp. 177-201. 
{ Cf. CUNNINGHAM, Proceedings of the London Mathematical Society, ser. 2, 
vol. 4 (1906), pp. 374-383. 
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Thus we reach a system of 2n linear equations of the first order in z to deter- 
mine the 2n functions 


Mis Mas Mais Mor °° 


This system of equations is analogous to (35). It has clearly singular points 
at 2=a,/2,z=4,/2,---,2=a,/2 and z= o, the last of which is a regular 
singular point. The finite singular points are not in general regular, although a 
slight consideration proves them to be of rank g = 0 at most. 

It is possible however to make a given one of these singular points z = a,/2 
a regular point by first making in (36) the simple preliminary transformation 


i i 


The new system of equations (36), and hence also (39) will be the same as 
before except that p‘') is replaced by p'')— 8 (i=1,2,---,n). These n sys- 
tems of equations (36), in general distinct and obtained by taking r=1,2,---,n, 
replace the single system (35) of the case q=9. 

The rth system of equations (39) has a regular singular point at z = a,/2. 
For the system (39) may be written 


dn, —1 
de 22 — + = + py (i=1, 2,---,m), 
dn, 1 n n n 
& “ia [2+ 2 t Piet) 
i = = 
A sufficient condition that z = a,/2 be regular is that there is no pole at z = a,/2 


of the second or higher order in the coefficients of 7,,, -+-, 7,45 Ms °**s Mao 
in the right-hand members. Clearly such a pole cannot exist except in the 
second n of these equations when i=7, and then only if p),+0. But 
in (39) p‘’ is replaced by p‘') — 8, which vanishes by (37). Hence the state- 
ment is proved. 

The exponents of the regular singular point z = a,/2 can be verified directly 


to be all 0 except two, which are 


Hence there exists a solution 7’), ---, ’), n°, ---, \ of the form 


id 2 id 2 ’ il 2 il 2 ’ 


where — a,/2) and 9)(z—a,/2) are analytic at z=a,/2. Each set 
of formulas 


—1-2 
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22 a; Nei Jj Nel a; 
1 


(i= 72, 


will then represent a solution provided that c, is a loop circuit to z= a,/2 so 
chosen that R[ x*(z— a,/2)] is negative along the ray. 

If (—/, —1)/2 is an integer, we can obtain solutions corresponding to the 
other exponent (— 7, — 2)/2 which will not be an integer. 

The proof of the validity of the above formulas is exactly like that in the 


r 


case g = 0. 
Let us now consider the general case when we have 


xq 


i ij 


1 
B,, (x) B, BY + 
We write 


fer" { Nig (2) + (2) (2) dz, 


where the »,,(z) are functions to be determined. First we have by substitution 


in (32) 
q n q q 
j=0 j=l k=1 t=1 


q 
+ > 2'n,, (z) dz (i=1,2,---,n). 
i=0 
The terms involve x’, x, ---,a*'. We propose to satisfy each of the equa- 
tions above by deleting in turn all but the terms in 
and aft!+* (s=0,1,---,q) 


and then satifying each of g +1 equations resulting. One obtains after the 
omission of a factor 2’ . 


J es 
j= +l=s 


k+l=q+8+ 
By an integration by parts 
fet att! w(2)dz = w(z) | — fe dw(z), 


in which the bracketed terms are again cancelled by a subsequent choice of 
limits, we may transform the coefficient e*'* in the integrands so that it is inde- 
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pendent of x. Equating these coefficients we then obtain in a formal manner 
the system of n(q + 1) differential equations for the functions 7,,, 


j=l k+l=7+8-+ j=l k+ 
(i=1, 2,---,n; s=0,1,---,q@). 


dyn, d 


These equations written in succession for s = q, q—1, ---, 0 are 
—[(¢q+1)z her Nig + (i 1,2,---,m), 


— ((¢+1)z—2,] - = de +t dz = 2 
j= 
(40) 


ij 


The singular points of (40) are at once seen to be z= a,/(¢ +1), z=4,/(q¢ +1), 

--,2=4,/(q¢+1) and z= o, the last of which is clearly a regular singular 
point while the others are of rank g at most. 

As in the case g=1 we can, however, make a given one of the points 

z=a,/(q¢+1) a regular singular point by making a preliminary transformation 


q 
The new system (32), and hence also (40), will not be changed except that 


Pi, pe”, p? ({=1, 2,-°+, 
will be replaced by 
py —8£,, per 1) - 

respectively. 

These n systems of equations (40) (r=1,2,---, 2) will, as in the cases 
q = 0 and q=1, lead to Laplace integral formulas. 

We shall now prove that the system (40) has a regular singular point at 
z=a,/(g+1). To prove this fact we note that, in view of the preliminary trans- 
formation assumed to have been made, the system (32) admits a formal solution 


T), To, ---, where 


TY =e Itly B.(#) (é= 3, 2, 


n 
Pall 
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Now it is easily proved that if the system (32) has a formal solution of this 
kind, then the adjoint system of equations, 


will have a formal solution 


1 


in which not every C, (j=1,2,---,) is zero. If this formal solution be 
substituted in the adjoint system, we obtain the following set of linear equations: 


=a,C; (j =1,2,---,), 


aC =a, Cpe (j=1,2,--+,”), 


From the first n of the above equations we infer that 
C,=0 (j+7), C_+0. 


Multiply now the nm equations of the ( 4 + 1 )th line in (40) by C,(i = 1, 2, ---, ») 
respectively, the set by (i =1, 2,---,), and so on. Then add the 
first nk equations thus obtained. In view of the above relations between the 
BY? one finds 


= 


| 


dz 


in which L,(7) denotes a certain expression linear in ;; with constant coeffi- 
cients. Since C,= 0 (i +7), C, + 0, we may write this in the form 

dy, 1 d 1 


r t=) 
(k= 0,1, 2,---,@). 


We are now able to prove that z = a,/(q +1) is a regular singular point. 
This will be proved if we show that when we solve (40) for the derivatives of 
n;; in terms of the ;; themselves, the coefficients of the expressions obtained 
have at z = a,/(q + 1) poles of the first order at highest. The expressions for 
dn,,[dz (i=1, 2, ---, m) certainly satisfy this condition, as the first set of 
equations (40) show. In solving for dn, _,/dz (i= 1,2, ---,) we distinguish 
the cases i+ r and i=r. Since in the first case dy, _,/dz appears with a 
coefficient —[(¢+1)z—a,] not zero at z = a,/(q+1), it is clear that poles of 


dz, 
(2), 
) 
n 


1909] ORDINARY LINEAR DIFFERENTIAL EQUATIONS 461 


higher than the first order do not appear. To solve for dn, ,__,/dz we use 
equation (41) for 4 = 1 and deduce the same fact. In this way we solve for 
dn, dn,,_,/dz, dn,/dz, each time making use of the appropriate 
equation (41), to prove that no poles of higher than the first order arise. Thus 
the statement is finally proved. 

The n(q + 1) exponents at z = a,/(q¢ +1) are all zero excepting ¢ + i of them, 
which are (— 1, — 1)/(q +1), (—4 — 2)/(q + + 
I have not found a direct proof of this fact, such as no doubt exists. An out- 
line of an indirect proof is found in the foot-note p. 462. 

We will have then a solution n''? (i=1,2,---,n; 7=0,1,---, q) such 


that 
a —/,-1 a 


where J) [z—a,/(q¢+1)] is a function analytic at z=a,/(g+1). This 
leads at once to the formulas 


where c¢, is a loop circuit about z=a,/(¢ +1) so chosen that 
R {at [2 —a;/(q¢+1)]} is negative along the ray. These actually repre- 
sent solutions, as can be proved just as in the cases g = 0 and g= 1. 

If (— 72; —1)/(q +1) is an integer, one may use solutions corresponding to 
the remaining non-zero exponents, none of which are integers. 


§ 5. Asymptotic representation. 
If one makes the substitution 


in the formulas (42) for 7 = r, there is obtained the expression 


—l,—1 


r 


where c; is a loop circuit to w=0. When one expands J,,(a~?-'w) and 
evaluates the successive terms by means of the formula 


there results for YY) a sum of terms of the same form as in the expression 7'Y) 
plus terms which approach zero as x becomes infinite.* 


* For a more complete consideration of a like point see the paper by Horn previously referred 
to, Mathematische Annalen, vol. 50 (1897), pp. 527-531. 
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In order to state the result concisely we shall modify slightly the conventional 
definition of asymptotic representation. 

Derinition. We may say that Y, is asymptotically represented by T’ 
along the ray arg x = >, or in symbols 


™, = ¢, 
if for some sector including arg «= ¢$ as an interior ray and for all values of m, 


AW ¢ 
AL + — +--+ + 


x |? 


in which lim e = 0 as x becomes infinite. 

The result is then that we have 

YP~ TP 

along any ray for which 2 { x**'[z—a,/(q+1)]} is negative. These expres- 
sions 7’/) are, of course, identical with the Z7’'Y’ which occur in the formal solu- 
tions of (32).* 

The ray of the loop circuit to z= a,/(q + 1) may vary in direction as long 
as it does not intersect one of the remaining x —1 finite singular points 
z=a,/(¢+1). It isa simple matter to give the actual sectors on which in 


general the above integral formulas hold. Denote by 7,, 7, ---, Ty the 
N=n(n—1)(q +1) rays in angular order for which 
R {(a,—a,)a*"} =0 (j+#k), 


and also write Ty,,; = 7, +27. We assume these rays to be distinct. Further- 
more, the real part of some particular one of the differences (a, — a,)a‘*! 
changes from positive to negative as we pass from the sector (7,_,, 7,,) to 
(7,5 Tns;)+ Denote the subseripts in this difference by j, and &,. The ray 
m m+1 m m 
of the loop circuit c; has an argument intermediate between a certain successive 
pair of arguments 
arg (%,,— a _,) and arg (a, — a, ) 


for which k,, =k, =j. In this case PR { x**'[z—a,/(q¢ +1)]} remains nega- 


tive for 
< ty +— 
The formula (42) 
VY, = PP 

* At this point it is easy to see that the exponents at z—a,/(q-+- 1) have the values above 
stated. For by comparison of the Laplace integrals and the formal solutions we get these values as 
possible exponents, and no others. But in the special case for which P,(z) =0, Py=p{; + az", 
it may be at once proved that the exponents have precisely these values. Then if the coefficients 
of the polynomials are changed continuously, the exponents are also so changed and hence retain 
the form specified. 
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represents a set of ¢g +1 solutions fixed by assigning the sector in which the 
ray of the loop circuit is to lie. This formula is valid for any one of the qg+1 
corresponding sectors 
+ q +1 ° 

It is to be noted that m and m’ are any pair of successive values of m for which 
ki =k, =j. 

If now we pass back to the solutions of (1), we infer that there exist solu- 
tiONS Y,5 Y, Of (1) such that 


y,~ SY (i=1, 2, ---, 
within any given sector 


< argx<tT., 
m™ 9 + q + 1 
This result has also been derived in certain special cases by Horn,* who has 
used the method of successive approximation. 


§6. Characterization at z=. 


At this point we can readily give a characterization of the nature of the solu- 
tions of (82) and thence of the system (1) equivalent to it at x = oo. 
THEOREM. There exist N=n(n—1)(q+1) fundamental systems of 
solutions of the canonical system (32), 
such that 
TY, Tn = arg < Tm+19 


and such that 
= VR (i+ 


Yon = + A, Pm, 


i,m+1 im 
provided we write = 


Proof. In order to avoid confusion we will state at the outset that the solu- 
tions (42) are not in general the solutions referred to in this theorem. The 
latter retain their asymptotic form throughout minimum sectors, while the former 
retain it throughout maximum sectors. 


*Crelle’s Journal, vol. 133 (1907), pp. 19-67 
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An immediate conclusion of the results of the preceding paragraph is that 
there exists for every ray arg a set of solutions 
(j=1,2,---,) such that along this ray 


ry ‘alias 


In virtue of this fact it is possible to separate the x-plane into some finite num- 
ber of closed sectors ¢, in each of which there exists a fundamental system of 
solutions V\,---, VP (j=1, 2,---,n) such that PY~ for x 
in o. We choose these sectors so that the rays arg «=7,, are included at 
most one within each sector c. 

Throughout a sector o not containing an inner ray 7, every solution 
Y,, ¥,, ---, ¥,, of (32) has a definite asymptotic form. For the general solu- 
tion of (32) is 


which gives the asymptotic relationship 
Ty. 


The relative magnitudes of these » terms as 2 becomes infinite in o clearly 
depend on the relative magnitude of 


(45) R(a,xt*'), ---, R(a at). 


This does not change except at the rays arg x= 7,,7,,---, Ty. Suppose for a 
moment that indices 1 to n are so chosen that the quantities (45) appear in order 
of descending magnitude. If ¢, is the first of the quantities c,, c,, ---, ¢, not 
zero, one has clearly for x in 

Furthermore, since the successive regions o have a ray in common, every 
solution Y,, Y,,---, ¥, preserves its asymptotic form throughout successive 
sectors o until one which contains a ray T,, 7,,---, Ty is reached. We may 
therefore unite into a single sector o, the region o which contain a given ray 7, , 
and all those which succeed it in angular order up to the one which includes 
the ray 7,,,, and affirm that on each of these sectors o there exists a 
fundamental system VY, VY), ---, PY) (j=1, 2,---,n), for which YY ~ T? 
throughout the sector. 

Next let us determine the character of the general solution VY, Y,,---, ¥,, 
in one of these sectors ¢,, which includes within it the ray arga2—=7,. Here 
again the relative magnitude of the quantities (45) determines the important 
term in (44). But R{P, (x)} and R{P, (x)} change order of magnitude 


along argx=T,, and moreover 
+1 
> 0 (argz<tm), 


{( a, — a, <0 (argz>Tm). 


1909] ORDINARY LINEAR DIFFERENTIAL EQUATIONS 465 


If we suppose indices so chosen that the quantities (45) are in order of magnitude 
at the first boundary ray of ¢, , two successive terms of (45), say the sth and the 
(s + 1)th, will change order of magnitude for arg2=7,,. We see by (44) 
that the solution Y,, Y,, .--, 7, will preserve its asymptotic form in ¢,, unless 


mc 


In this case the solution is of the asymptotic form 
at the initial ray of the sector, and of the asymptotic form 


at the terminal ray. 


If, however, we have two solutions such that at the initial ray of the sector 


stl 


(¢=1, 


we can always choose a combination Y = Y,+ AY; so that (¥\’, 


preserves its asymptotic form c,7'\’, c,7S’,---,c,Z throughout the sector 
co. For we have 


m 
7 (s+1) (8) 


Y,=c' 


i stl 


Fy +e 


The required value of A is —c,,,/c’ 


‘at1° 


Let now V(), VSP, ---, (j=1, 2, ---, m) be any fundamental sys- 


ni 


tem of solutions of (32) such that 


for x on o,. All these solutions will preserve their asymptotic form throughout 
o, except perhaps the solution 


Pips 


in view of the preceding discussion. However, in this case we can choose a con- 
stant A, so that the solution 


rts k j 7 (k (k 
preserves its asymptotic form 
TH 


throughout the sector ¢,. Therefore if we write 


YH YH A,Y* 
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the new fundamental system of solutions ---, (j=1, 2,---, 2) 
preserves its asymptotic form throughout the sector o,. 

Passing now to the sector o,, we obtain a third fundamental system of solu- 
tions, and so on. Thus we define ---, successively in terms 
of I‘). 

These fundamental systems of solutions have therefore the desired properties 
stated in the theorem except that we need not have 


These relations, too, will be satisfied when a proper choice of the initial funda- 
mental system 
(46) YY, ---, FY (j=1, 2,---, 
is made. To prove this fact we show first that the final fundamental system 
of solutions 
(47) PP ysis PP (fj =1, 2,---, 2) 
is entirely independent of the choice of the initial system (46). 
Let 
Ze, Zi, (j= I, 2, m=1, 2, N) 


be a second set of solutions of the same kind; the constants which correspond 
to A, will be denoted by B,,. 

Let R(a,x**') be that one of the quantities (45) which is least in 
magnitude in(7,,7,). There can only be one solution of (82) asymptotic to 
Ty, +--+, T on (7,, 7,), and we must have 


(48) YY? Zi (i=1, 2,---, 2) 


in this sector. Now each pair of the quantities (45) become equal for 2(g + 1) 


values of 7, 


It is clear, then, that on the sector 


each pair of these quantities (45) become equal once and only once. The quan- 
tity W(a,2**') will therefore increase over the remaining quantities (45) as 
arg « changes from 7, to 7,-+ 7/(q¢+1). It is therefore clear that 

= 


iu? 


ZY = ZY ZY. 


T —1 
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Hence we find 


for all values of m less than uw + 1. 
Next let R(a,x**') be that one of the quantities (45) which is second in 
magnitude in (7,,7,). We must have in this case 


7 Y(t) 
Z = I + el} 


for m= 1, 2, ---, until a value m= @ is reached for which /2?(a,x**') increases 
over R(a,x?*'). It must be remembered that /?2(a,x”**') increases over all 
the quantities (45) except R(a,x**') in (7,, 7,). Now in virtue of the 
relations 


(t t) (8) *) *) 
= 234+ B22, = 


i, 6+1 


the preceding equality for m = 6 may be written 
Also since 2(a,x‘*') is greater than 2?(a,x**') in (7,, T),,) we conclude that 


—B,=—A,+¢. 
We have already shown that 
and we obtain from (50) 
Since for m > 6, the quantity (a,x**') will increase over the remaining 
quantities (45), we infer that 
YQ, = 


— @ 
im im (i= 1, 2, n) 


for all values of m greater than @. 
Likewise we have 
where 2(a,x*t') is the third of the quantities (45) in order of magnitude. 
The relation 


im im 


will hold for successive values of m, the constants c and d only changing for a 
value of m such that one of the three quantities 


R(a,x**"), R(a,x%*"), R(a,x"), 


changes order of magnitude at argz=7,. If the first increases over the 
second, the constant d may change; when the second increases over the third, 
the constant c becomes zero; when the first increases over the third, the con- 


ry ZS (i= i. 2, n) 
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stant c becomes zero. From and after this fixed value of m < wu + 1, for which 
the first increases over the third, we have 


PQ) = (¢=1,2,---,m). 


In this way we prove that after some fixed value of m < + 1 we shall always 
have 
=_ ZY) ({=1, 2,-->, n), 


tm am 


whatever be the particular value of 7. Therefore we have 


im 


VD ZY (4, j=1, 2, #) 


for m = pw, and of course for greater values of m, in particular m = V+ 1. 

Choose Y'‘j) to be this fixed set Y\/}.., multiplied by e~**4*—", and all the 
conditions of the theorem will be satisfied. This is a possible choice, since 7\ 
after a complete circuit of «2 = co becomes ¢?*4*-!7'), and therefore we have 
~ in (7,,7,). This completes the demonstration. 

The theorem whose demonstration has just been completed leads at once to a 
like theorem concerning the solutions of (1). 

TueoremM. There exist NV = n(n —1)(q + 1) fundamental solutions of (1), 


(1) (1) 


(1) 
Yim? Yom? Fan? 


(n) (n) (n) 


Yim? Yom? 
such that 


oD) 
and such that 
Ye = Yim (i+ jm), 
(Kin 


provided we write = 
This is a complete characterization of the solutions of (1) in the vicinity of 
x = 0, since any system of functions ¥’'’) satisfying these conditions will form 


the solutions of a system of equations (1). The characteristic constants which 
determine the nature of the solutions are the n(g + 2) exponential constants 


Ay, K, (j= 1,3 n), 


and the n(n —1)(q +1) transformation constants A,, A,,---, Ay. Not all 
of the constants A, are independent, since instead of the particular n formal 
solutions used we might have multiplied each of these by an arbitrary constant. 
The new transformation constant would then be altered to 


Cy 
A,,= c 
Jim 
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Two sets of transformation constants thus related are considered to be the same. 
Any two systems (1) and (1) with the same characteristic constants (except that 
x, and %, may differ by integers) may be proved to be equivalent at x = 00, and 
conversely. 

The total number of essential characteristic constants is accordingly 


But we can dispose of n — 1 further constants of (82) by the transformation 


(¢=1, 2, ---, 
It is clear then that the canonical system of equations (82) also has at most 
n?(q +1) +1 essential constants. If it be admitted that no further reduction 
of (32) is possible, we infer that the characteristic constants are not connected 
by any necessary relation. 


$7. Generalized Riemann problem. 


We now propose the following generalization of the Riemann problem: To 
construct a system of n linear differential equations of the first order with pre- 
scribed singular points 

of respective rank 

and with a given monodromic group, the characteristic constants being assigned 
for each singular point. 


Consider the system of equations 


m U+2 4 


j=1 \k=1 \4 


Ym+1 
ijkl —_ 


with singular points x,, x,,---, = 00 of prescribed rank ¢,, 9.5 


m? ““m+1 


In this equation appear 


2m + 1) 


f=, 


arbitrary constants. 

In order to fix the problem we consider the fundamental system of solutions 
to 6,,, and take the group of this particular fundamental system of solutions 
to be assigned. 

The monodromic group contains n*m arbitrary constants, n* for each funda- 
mental substitution. Furthermore, the characteristic constants are 


41] 
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in number; but between these constants and those of the group there are 
n(m +1) relations, since the n exponents at each of the m + 1 singular points 
are determined either by the group or by the characteristic constants. 

We must further prescribe at each singular point precisely which solutions 
are to correspond to the canonical fundamental system of solutions given in the 
second theorem of § 6. Now the group itself determines the correspondence 
(since it determines the exponents) except for n multiplicative constants. There 
are accordingly 


(n—1)(m+1) 


further conditions. Thus there are in all 


>, 
conditions to be satisfied. 
The number of conditions equal the number of constants at our disposal. 
This is the basis on which we postulate the possibility of the general problem. 
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AUTOMORPHISMS OF ORDER TWO* 


BY 
G. A. MILLER 


§ 1. Introduction. 


If an operator ¢ transforms a group G@ into itself without being commutative 
with each operator of G, and if ¢* is commutative with each operator of G’, then 
t is said to transform the operators of G according to an automorphism of order 
two. Moreover, every automorphism of G which is of order 2 may be obtained 
by transforming G' by operators having the given properties of ¢.¢ Let s be 
any operator of G and assume that 


t—'s,t=8,8,. 


Since ¢? is commutative with every operator of G it results that 


—lg gf = as 
t's, st =8,8;s=8, or 8,8, = = t's, t. 


That is, in every automorphism of order 2 some operator besides the identity 
must correspond to its inverse. It also results that the necessary and sufficient 
condition that an automorphism shall be of order 2 is that every commu- 
tator arising from the automorphism corresponds to its inverse under the 
automorphism. 

It is well known that automorphisms play a fundamental role in many prob- 
lems in group theory, but little has been done toward a general theory of auto- 
morphisms. Even the automorphisms of order 2 present many difficulties at 
the present stage of development of group theory. Such a fundamental question 
as whether every group admits automorphisms of order 2 has not yet been 
answered. When G is an abelian group whose order exceeds 2 it clearly 
admits at least one automorphism of order 2. 

The main results of the present paper relate to the case when G is abelian. 
In § 4 a few results are obtained which relate also to non-abelian groups, and 
the final theorem is an extension of results coming under the heading of the 
present paper. 


* Presented to the Society (Chicago), April 9, 1909. 
t FrRopentius, Berliner Sitzungsberichte, 1895, p. 185. 
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§ 2. Abelian groups of odd order. 


When @ is an abelian group of odd order the identity is the only commu- 

tator arising from ¢ and invariant under ¢. All the commutators which arise 
from transforming G by ¢ constitute a group 7, and to each operator of the 
quotient group G/H = H’ there corresponds one and only one operator of G 
which is commutative with ¢. The independent generators of G may be 
selected from /7 and these invariant operators under ¢. Hence the theorem: 
_ The independent generators of an abelian group of odd order may be so 
chosen that in any given automorphism of order 2 each of these generators 
corresponds either to itself or to its inverse. In particular, if an operator of 
order 2 transforms into itself an abelian group of odd order, it is possible to 
select the independent generators of this group in such a manner that each of 
them is transformed either into itself or into its inverse. 

If a set of independent generators of G is so selected that the order of each 
of them is a power of a prime, every other possible set of such independent gen- 
erators may be placed in a (1, 1) correspondence with this set. Hence it is 
easy to find the number of sets of conjugate operators of order 2 in the group 
of isomorphisms (J) of G. This may be accomplished as follows: Represent 
all the independent generators, which have been selected in the given manner, 
by letters in such a way that the generators of the same order are represented 
by the same letter and the generators of different orders by different letters. 
Form all the possible combinations of these letters, taking one, two, three, ete., 
at atime. The number of these combinations is the number of distinct sets of 
conjugate operators of order 2 in 7. In particular, when G is of order p* and 
of type (1, 1,1, ---) the number of complete sets of conjugate operators of 
order 2 in J is a. 

By means of the results obtained in the preceding paragraph it is easy to 
determine the number of possible groups of order 2(2n + 1) involving a given 
abelian subgroup of order 2n +1. In fact, the number of these groups is 
exactly one more than the number of*complete sets of conjugate operators of 
order 2 in the J of the abelian group of order 2x + 1. In particular, there are 
exactly a + 1 groups of order 2p*” involving the abelian group of order p*” and 
of type (m, m,m,---). Since the group of isomorphisms of a Sylow subgroup 
of G@ involves just one invariant operator of order 2,* the total number of 
invariant operators of order 2 in I is 2*—1, X being equal to the number of 
distinct Sylow subgroups in G'. 

For the purpose of determining the total number of operators of order 2 in 
I it is convenient to observe that the group of isomorphisms of any abelian 
group is the direct product of the groups of isomorphisms of its Sylow sub- 


*Transactions of the American Mathematical Society, vol. 2 (1901), p. 260. 
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groups. If m,, m,, m,,--- represent the numbers of operators of order 2 in 
the groups of isomorphisms of the various Sylow subgroups of G', the number 
of operators of order 2 in J is (m, + 1)(m,+1)(m,+1)---—1. Hence we 
may confine our attention to the determination of the number of operators of 
order 2 in the group of isomorphisms of an abelian group (/) of order p”, p 
being an odd prime. When P is cyclic its group of isomorphisms involves only 
one operator of order 2, viz., the invariant operator of this order, corresponding 
to the isomorphism in which every operator of P corresponds to its inverse. 
In all other cases this group of isomorphisms involves non-invariant operators 
of order 2 and the number of distinct operators in a complete set of conjugates 
is evidently equal to the number of conjugates under J of the corresponding 
pair of subgroups /7, /’. 

In general, 7, HT’ are any two independent subgroups of G whose product 
equals G. Hence the number of operators of order 2 in J is equal to the num- 
ber of pairs of independent generating subgroups contained in G, when /7, H’ 
and #7’, H are regarded as distinct pairs, and #7 differs from the identity 
while #7’ may be the identity. This number is clearly equal to the sum of the 
quotients obtained by dividing successively the order of J by the products of the 
orders of the groups of isomorphisms of /7, 1’ when 7, H’ represent succes- 
sively all the possible types of pairs of independent subgroups (except that // is 
not the identity) contained in G. When // is the direct product of Sylow sub- 
groups, the corresponding quotient is equal to unity and vice versa. Hence this 
is another way of seeing that the number of invariant operators of order 2 in 
T is equal to 2‘ — 1, A being the number of Sylow subgroups in G. Some of 
these results may be expressed in a convenient form as follows: The number of 
operators of order 2 in the group of isomorphisms (I) of an abelian group 
( @) of odd order is equal to the sum of the quotients obtained by dividing 
successively the order of I by the products of the orders of the groups of iso- 
morphisms of all the possible pairs of independent generating subgroups H, 
HT’ of G such that H is not of the same type in two pairs and that H is not 
the identity. Each of these separate quotients represents the number of oper- 
ators of order 2 in a complete set of conjugates under I. 

To illustrate this theorem we shall employ it to determine the number of 
operators of order 2 in the well known group of isomorphisms of the abelian 
group of order p*(p > 2) and of type(1,1,1). The order of J is evidently 
Pp (p-1 )(p? —1)(p—1), and H, H’ are of one of the following three 
types (1), (1,1);(1,1),(1);(1,1,1),(0). The number of operators of 
order 2 in J is therefore equal to the sum of the following three quotients : 


P(p?+pt1), t+ptt), 1. 


Hence this J involves two sets of p’(p* + p + 1) conjugate operators of order 
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2, and, in addition to this, it involves the operator of order 2 which transforms 
every operator of G into its inverse. It may be added that H, H’ are neces- 
sarily products of Sylow groups whenever G is cyclic. This establishes contact 
between the given theorem and the well known result that the group of isomor- 
phisms of a cyclic group is abelian. As the type of H’ is completely deter- 
mined by that of H7, we may confine our attention to the latter in determining 
the number of sets of conjugate operators of order 2 under J. 


$3. Abelian groups of even order. 


When the order of G is the double of an odd number its group of isomor- 
phisms is the same as if its order were this odd number. This is a special case 
of the theorem that the group of isomorphisms of G is the direct product of 
the groups of isomorphisms of its Sylow subgroup. Hence we may confine our 
attention to the case when the order of G is of the form 2”. This case presents 
much greater difficulties than the one considered above. As every non-cyclic 
abelian group of order 2” is the direct product of abelian groups whose invari- 
ants are equal to each other, we shall generally confine our attention to the 
determination of automorphisms of order 2 in such an abelian group. Through- 
out the rest of this section it will be assumed therefore that the order of G is 
2” and that all the invariants of G are equal unless the contrary is stated. 

Representing the commutator subgroup of G by /, as was done in the pre- 
ceding section, it is easy to prove that each of the independent generators of /7 
must be of order 2, 2*~', or 2*, when 2* (a> 1) is the common order of the 
independent generators of G. In fact, since each operator of H/ is transformed 
into its inverse by ¢ and since each operator of G is generated by an operator 
of order 2* contained in G‘, an operator (s) of 7, whose order exceeds 2 , must 
be generated by an operator of order 2* which generates only two operators that 
are commutative with ¢. This generator of s is transformed by ¢ into itself 
multiplied. by an operator of order 2*~', since ¢—'st = s~*s. As this operator 
of order 2°~' generates s*, it results that the square of every operator of H 
whose order exceeds 2 is generated by an operator of order 2*-' which is con- 
tained in H. This proves that the independent generators of H whose orders 
are divisible by 2°—' generate the squares of all the operators of H whose orders 
exceed 2. Hence there are three cases to consider: 1) When all the operators 
of H except the identity are of order 2, 2) When // involves operators of 
order 2*~' but none of order 2*. 3) When / involves operators of order 2. 
These three cases will be considered in the given order. 

When // involves no operator of order 4, the order of H must divide 2’, 
d=m/a. The number of automorphisms of order 2, when the order of H is 
equal to 2%, is equal to the order of the group of isomorphisms of 47. In gen- 


— 
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eral, when the order of H is 2*-', 7 <d, the number of these automorphisms is 
equal to the order of the group of isomorphisms of // multiplied by the product 
of the number of the subgroups of order 2' in the group generated by all the 
operators of order 2 in G and the number of subgroups of order 2“ and of type 
(a, a,a,---) contained in G. As the numbers of these subgroups are well 
known,* the problem of determining all the automorphisms of order 2 that arise 
in the given manner is completely solved. Hence the theorem: Jf G is an 
abelian group of order 2"* and of type (a, a, a,---),#>1, the number of 
automorphisms of order 2 which give rise only to commutators of order 2 
is equal to the sum of the products obtained by multiplying the order of the 
group of isomorphisms of the abelian group of order 2*~' and of type 
(1,1, 1,---) by the product of the number of the subgroups of order 2' and 
of type (1, 1, 1, ---), and the number of subgroups of order 2" and type 
(a, a, a, ---) contained in G, for every value of l= 0,1, ---,d—1. 

When a=1, both G and H are of type (1, 1, 1, ---) and hence each 
operator of 7 is commutative with ¢. From this it follows that the number of 
independent generators of H is equal to or less than half the number of these 
generators in G. It is also evident that the number of sets of conjugate 
operators of order 2 in J is equal to the number of possible orders of /7; i. e., 
it is equal to the largest integer which does not exceed d/2. If the order of H 
is 2 the number of operators of /7 which are commutative with ¢ is 2“~‘ and 
the number of operators of J which are conjugate with an operator correspond- 
ing to such an // is clearly equal to the product of the following three numbers: 
1) The number of subgroups of order 2“~* in G. 2) The number of sub- 
groups of order 2’ in such a subgroup of order 2“~*. 3) The order of the 
group of isomorphisms of #7. The total number of operators of order 2 in 
the group of isomorphisms of the abelian group of order 2“ and of type 
(1, 1, 1, ---) may therefore be found by addition of these products for every 
value of A from 1 to the largest integer which does not exceed d/2.¢ For 
instance, the group of isomorphisms of the abelian group of order 2‘ and of type 
(1, 1, 1, 1) contains two sets of conjugate operators of order 2 involving respec- 
tively 105 and 210 operators, as may easily be verified, since this group of 
isomorphisms is the alternating group of degree 8. 

The simplest case in which 7 involves operators of order 2*-' is that in which 
every operator of G corresponds to its inverse. In this case it is evident that 
HT is the direct product of d cyclic groups of order 2*-! and that the corre- 
sponding operator of order 2 in J is invariant. When a= 2, this is included 
in the preceding two paragraphs. Hence we may assume a> 2. Since the 
given automorphism corresponds to an invariant operator of order 2 in J, it 


*Cf. Annals of Mathematics, ser. 2, vol. 6 (1904), p. 5. 
t Bulletin of the American Mathematical Society, vol. 8 (1902), p. 391. 
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may be followed by each one of the group of automorphisms giving rise to com- 
mutators of order 2 only. If it is followed by any other automorphism G 
must involve invariant operators of order 4 under the automorphism, since each 
commutator arising from the automorphism corresponds to its inverse in the 
automorphism. Hence the theorem: The number of the automorphisms of 
order 2 in which H involves all the operators of order 2*—', a> 2, contained 
in G is one more than the number of automorphisms of H which are of order 
2 and lead only to commutators of order 2. As the latter number was con- 
sidered above, the determination of these automorphisms is complete. 

Suppose that /7 involves an operator (c,) of order 2*. We proceed to prove 
that G must contain an invariant operator of order 2* in this case. In fact, if 
t-'s,t =c,8, it is necessary that s, be of order 2* and also independent of ¢,. 


Hence the group generated by c, and s, isof order 2°. Since ¢~'c,t=c;'=c7*e 


and ¢~'s?t = c? s?, it results that c,s; is invariant under ¢, and it is evident that 
c,s; is of order 2*. This proves the theorem: When // involves an operator 
of order 2* an operator of order 2* in G must correspond to itself in the auto- 
morphism of order 2 which gives rise to H. It is also evident that c, s? and ¢, 
generate groups which have two and only two common operators. 

The group of order 2** { ¢,, s, } cannot involve more than one commutator of 
order 2 since these commutators are invariant under ¢. It results from this 
that if H involves two independent generators {c,, c,} of order 2°, the group 
generated by c, has only the identity in common with {c,, s}, as the latter 
involves no commutators under ¢ except such as are generated by c,. Hence ¢ 
must transform a fourth independent generator s, as follows: ¢~'s,t = c,8,, 
whenever /7 involves at least two independent generators of order 2*. This 
method of argument may be repeated until all the independent generators of 
order 2* in #7 have been exhausted. As these arguments apply to any abelian 
group of even order the result may be stated as follows: Jf the commutator 
subgroup arising from an automorphism of order 2 of any abelian group of 
even order ( & ) involves n independent generators of highest order in K, then 
K involves at least 2n independent generators of this order. 

It was observed in §2 that every abelian group of odd order has an auto- 
morphism of order 2 giving rise to a commutator subgroup which involves all the 
independent generators of the abelian group, while the preceding theorem shows 
that the commutator subgroup arising from an automorphism of order 2 of an 
abelian group of even order cannot involve more than one half of the largest 
independent generators of the group. If this commutator subgroup involves 
exactly half of the independent generators of highest order, it cannot involve 
more than half of those of the second highest order whenever this order is also 
even, ete. All the operators of which are not powers of operators of higher 


order contained in /7 must arise from operators of highest order in G. 
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§ 4. Commutators of order p*. 


In the preceding section it was observed that an automorphism of order 2 of 
an abelian group of even order must always give rise to at least one invariant 
commutator of order 2. Moreover if all the commutators arising from a given 
automorphism are of order 2 and correspond to themselves in the automorphism, 
this automorphism must be of order 2. In the present section we inquire into 
the properties of a non-abelian group (A) which has no commutator whose 
order exceeds 2 but which is not restricted in any other way. It is not assumed 
in this section that all the commutators of A arise from automorphisms of 
order 2. 

Since all the commutators of A are of order 2, every operator (s) of order 2 
in A must be transformed under A either into itself or into cs, where c is of 
order 2. Moreover, from the fact that both s and cs are of order 2 it results 
that c, s, and cs are three commutative operators which, together with the 
identity, constitute the four-group. That is, every complete set of conjugate 
operators of order 2 in A’ is composed of mutually commutative operators. In 
other words, all the operators of order 2 in A are found in invariant subgroups 
each of which is of type (1,1,1,---). From this it follows directly that A 
involves only one Sylow subgroup of order 2”, and that the quotient group with 
respect to this Sylow subgroup is abelian. This result will be generalized in 
the following paragraph. 

If the orders of all the commutators of a non-abelian group K are of the form 
ps, p, being a prime number, it follows directly that every Sylow subgroup of 
K whose order is not a power of p, is abelian. Moreover every operator of A’ 
that transforms into itself a subgroup of order pz, (a +1), is commutative 
with every operator of this subgroup. Hence it results from a theorem due to 
FroBEntus that A contains a characteristic subgroup of index p°®, where p? is 
the highest power of p, that divides the order of A.* As the orders of the 
commutators of this characteristic subgroup must also be of the form p*, we may 
repeat the given argument when the order of this subgroup is not itself of the 
form p*. Hence the following theorem: Jf the orders of all the commutators 
of a non-abelian group K are of the form p>, p, being any prime number, K 
involves only one Sylow subgroup of order p and the corresponding quotient 
group is abelian. 

From the preceding paragraph it follows directly that every group in which 
the orders of all the commutators are powers of the same prime number must be 
solvable. It is also easy to see that A contains an abelian subgroup which 
is the direct product of Sylow subgroups of every order except p”’. This 
follows directly from the fact that all the Sylow subgroups of order p® are 


* FROBENIUS, Berliner Sitzungsberichte, 1901, p. 1324. 
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contained in a characteristic subgroup of order p’'p*. That such a group as 
K may involve more than one Sylow subgroup of every order except p’’, results 
directly from the fact that the metacyclic group of degree p and of order ° 
p(p—1) is included among those which satisfy the conditions imposed on X. 


URBANA, ILLINOIS. 


RESOLUTION INTO INVOLUTORY SUBSTITUTIONS OF THE TRANS- 


FORMATIONS OF A NON-SINGULAR BILINEAR 
FORM INTO ITSELF* 


BY 


DUNHAM JACKSON 


In a recent number of the Transactions of the Connecticut Acad- 
emy, E. B. Wixson obtains a necessary and sufficient condition that a linear 
transformation be factorable into two involutory transformations. If this condi- 
tion is translated from the language of dyadics, which WILSON uses, into that of 
the algebra of matrices, its form suggests at once a theorem of FROBENIUS con- 
cerning transformations of a bilinear form into itself. A mere combination of 
these two theorems is sufficient to establish the following, which is simpler in 
statement than either : 

I. A necessary and sufficient condition that a linear transformation be such 
as to carry some non-singular bilinear form into itself, is that it be factorable 
into two involutory transformations. 

It will be seen that the part of this theorem relating to the necessity of the 
condition is a generalization of a theorem of P. F. Smira + which WILg0N uses 
in proving his theorem, and which is substantially as follows : 

II. (SmitH). A necessary condition that a linear transformation be such 
as to carry some non-singular quadratic form into itself, is that it be factor- 
able into two involutory transformations. 

The converse of SMITH’s theorem is not true, as is remarked by WILSON; a 
substitution which is factorable into two involutory transformations is not neces- 
sarily capable of carrying a non-singular quadratic form into itself. In fact, 
the question of the possibility of factoring a transformation into two involutory 
transformations has no essential relation to the subject of symmetric bilinear 
forms, as distinguished from bilinear forms in general. $ 


* Presented to the Society (Princeton), September 13, 1909. 

+ P. F. Smita, these Transactions, vol. 6 (1905), p. 13. 

t For a necessary and sufficient condition that a transformation be capable of carrying a non- 
singular quadratic form into itself, see FROBENIUS, Crelle’s Journal, vol. 84 (1878), p. 41. 
Compare with this the theorem numbered IV below. 
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The theorems of W1Lson and FROBENIUS mentioned at the beginning of the 
article may be stated thus: 

III. (Wizson).* A necessary and sufficient condition that a linear trans- 
formation be factorable into two involutory transformations, is that each 
elementary divisor of its characteristic matrix be paired with another, whose 
degree is the same, and whose root is the reciprocal of the root of the first, 
except that an elementary divisor which is a power of (>%—1) or (X+1) may 
be paired with itself. 

IV. (Fropentus.)¢ A necessary and sufficient condition that a linear 
transformation be such as to carry some non-singular bilinear form into itself,t 
is that each elementary divisor of its characteristic matrix be paired with 
another, whose degree is the same, and whose root is the reciprocal of the root 
of the first, except that an elementary divisor which is a power of (AX —1) or 
(~ +1) may be paired with itself. 

Theorem I is an obvious consequence of III and IV. But while Fropentus 
in proving IV uses the method of the algebra of matrices, W1ILson’s proof of 
III depends on Grpss’s method of dyadies and double products, and also on the 
work of Smirx, who uses still another method. It is possible, however, to give 
a proof of WILson’s theorem depending only on the algebra of matrices, and 
so to obtain greater uniformity of treatment. For a part of the proof it is nec- 
essary only to rewrite WiLson’s dyadics in the matrix notation; the use of 
SMiTH’s theorem in establishing the sufficiency of the condition is avoided by 
actually writing down the factors in all, as WILSON does in some, of the special 
vases that are first considered. 

In the course of the demonstration, reference will be made to another theorem 
of Frosentvs, slightly more general than the one already given: 

V. (Fropentus.)§ A necessary and sufficient condition that two linear 
transformations a, B be such as to carry some non-singular bilinear form 
into itself,|| is that the elementary divisors of the characteristic matrix of B be 
obtainable from those of a by replacing each root of the characteristic equation 
by its reciprocal. . 

From this, which is itself not difficult to prove, Theorem IV follows immedi- 
ately. For if it is assumed that the two sets of variables in the bilinear form 
are to be subjected to the same transformation, then 8 = a’, where a denotes 
the conjugate of a, and the theorem just stated reduces to IV, since the ele- 
mentary divisors of the characteristic matrix of a’ are the same as those of a. 

*E. B. Wiuson, Theory of Double Products and Strains in Hyperspace, Transactions of the 
Connecticut Academy of Arts and Sciences, September, 1908; p. 41. 

Tt Loe. cit., p. 34. 

t That is, when applied to each of the two sets of cogredient variables. 


Z Loe. cit., p. 31. 
|| That is, as FROBENIUS uses the terms, that a non-singular matrix ¢ exists such that a¢3 —0. 
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To proceed with the proof of WILson’s theorem: Suppose at first, that a is 
the matrix of the product of two involutory transformations, so that 


a= oy where = I, Y= I. 
a= ($y) = 


If 8, y are any two non-singular matrices, the elementary divisors of the 
characteristic matrices of Ay and of yf are the same, since 


yB = 


(see, for example, BOcHER, Introduction to Higher Algebra, p. 286). 

If y is any non-singular matrix, the elementary divisors of the characteristic 
matrix of y~' can be obtained from those of y by replacing each root by its 
reciprocal, by Theorem V above, for 


Then 


since 


yyy" 


In the present case, a, 6, are necessarily non-singular, since ¢* = y* = J. 
Therefore the elementary divisors of the characteristic matrix of a~' 
one hand, the same as those of a, and, on the other, the same as those of « with 


are, on the 


each root replaced by its reciprocal ; and this amounts to saying that the condi- 
tion stated in Theorem III is necessary. 

Next, suppose that a is any matrix such that the elementary divisors of its 
characteristic matrix satisfy this condition. The form of the condition insures 
that « is non-singular. It is to be shown that a is the product of two square 
roots of the unit-matrix. It will be convenient to make use of the following 

Lemma. If ais a matrix expressible in the form dy, where ¢? = = J, then 
any matrix §&, whose characteristic matrix has the same elementary divisors as 
that of a, is so expressible. 

For 8 may be written in the form 


B= yay = = 


—1 


where 


Tr 
=z J 


= voy" = 1? = YY 


= J. 


The truth of the following statements becomes obvious on writing down the 


and similarly 


formule : 

If the n-rowed matrix 8 has elements different from zero only in an m-rowed 
principal sub-matrix 8, and the n-rowed matrix y has elements different from 
zero only in the corresponding m-rowed sub-matrix ¥, then the product Sy has 
elements different from zero only in the corresponding m-rowed principal sub- 
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matrix, and these elements are those of the matrix Ay obtained by multiplying 
together 8 and 7 as m-rowed matrices. 

If the n-rowed matrices 8, y have elements different from zero only in the 
principal sub-matrices 8, 7 (not necessarily of the same order), neither of which 
contains any row or column corresponding to a row or column of the other, then 
By = 0. 

Therefore, if B=8,+8,+---+8, where the 
non-vanishing elements of §,, ---, 8, form non-overlapping m,, ---, m,-rowed 
principal sub-matrices, and those of ¥,, ---, y, occur in the corresponding sub- 


matrices, then 
By BY + + + Bo,» 


where in forming the products 8,y,, as far‘as the non-vanishing elements are 
concerned, the factors may be regarded as matrices either of order m, or of 
order n.* 

In particular, if 8 = 8, + 8,+---+ B, as before, and the sub-matrix in each 
8, is factorable into two square roots ¢;y, of the unit-matrix of order m,, then 
8 is factorable into two square roots of the unit-matrix of order 7, for 


It will be shown next that the condition of the theorem is sufficient in the 
special cases that the characteristic matrix of the given matrix has just one ele- 
mentary divisor (A+ 1)", or just two elementary divisors (A — a)’, (A—1/a)’. 
By the lemma above, it will be enough to exhibit in each case a single matrix 
which can be factored into two square roots of the unit-matrix, and whose char- 
acteristic matrix has the elementary divisor or divisors in question. The work 
that has just been done will then make it easy to complete the proof of the the- 


and 


orem in general. 

Let two matrices ¢, y be defined as follows: The elements in the principal 
diagonal of each shall be alternately 1 ard — 1, and those just above the prin- 
cipal diagonal alternately 1 and 0, the 1’s of ¢, above the diagonal, correspond- 
to the 0’s of y; all the other elements are zero. If n, the order of the matrices, 
is 4, to take a definite case, the formule are 


1 1 0 0 /1 0 0 0 

0-1 0 0 lo -1 1 

lo 0 O- 


* FROBENIUS, loc. cit., p. 18. 
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On multiplication, it appears that ¢? = J, y’ = IJ, and ¢y is a matrix whose 
characteristic determinant is equal to (> — 1)"; when n = 4, as above, 


1-1 1 0 
0 0 
0 0 1 -1 
0 0 0 1 


The first minor in the upper right-hand corner of the characteristic determinant 
does not vanish for X =1, so that there is just one elementary divisor (> — 1)". 
Therefore any matrix whose characteristic matrix has only a single elementary 
divisor, of the form (A — 1)", is factorable into two square roots of the unit- 
matrix. A similar result may be deduced for the case that there is a single 
elementary divisor (A + 1)", by considering the product (— ¢)¥. 
Let the notation 
a 
y 


represent the matrix of order 2e having the e-rowed matrices a, 8, 7, &, in its 
corners, arranged as indicated. A similar notation might be employed for 
matrices of order ne. By actually following through the process of multipli- 
cation, it is fairly easy to see that the ordinary rule for matrix multiplication 
still holds, if the elements of the matrices are regarded, not as scalars, but as 
themselves representing matrices, all of the same order.* 

Let $, + be two square roots of the unit-matrix of order e, whose product 
dy is such that its characteristic matrix has just one elementary divisor, (A — 1)’. 
It has been shown that such matrices exist. Let a be any number different 
from zero. Let 

0 V 0 
?= 1. ’ y= 


Vay 0 


where a may denote either square root, provided it is given the same value 
throughout. Form the product of ¢ and y: 


* Regard each matrix as a sum of matrices, in each of which all the matrix-elements but one 


are zero. 


| 
| 
oy = 1_.||- 
0 
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The characteristic matrix of apy has just one elementary divisor (A — a)’, and 
that of 1/a¢y has just one elementary divisor (A —1/a)’. Therefore the 
characteristic matrix of ¢y has just two elementary divisors, (A — a)’ and 
(A—1/a)y. Furthermore, 


where J represents the unit-matrix of order e, and J, that of order 2e. Con- 
sequently any matrix whose characteristic matrix has just two elementary 
divisors (A — a)’ and (Aj —1/a)’ is factorable into two square roots of the 
unit-matrix. 

Return now to the matrix a, which was assumed to be such that the elemen- 
tary divisors of its characteristic matrix satisfy the conditions of the theorem. 
It is possible to write down a matrix 8, having these same elementary divisors 
in its characteristic matrix, and expressible in the form 8, + 8,+ ---+8,, 
where the non-vanishing elements of £,, ---, B, form non-overlapping principal 
minors, and the characteristic matrix of the sub-matrix corresponding to each 8, 
has either just two elementary divisors, (A — a)" and (A — 1/a)*, or just one 
elementary divisor, (A-+1)". In either case, the sub-matrix is factorable into 
two square roots of the corresponding unit-matrix. Therefore 8, regarded as a 
linear substitution, is the product of two involutory transformations, and hence 
the same is true of a. 
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ON SINGULAR POINTS IN THE APPROXIMATE DEVELOPMENT 


OF THE PERTURBATIVE FUNCTION 
BY 
FRANK WALKER REED 


$1. Inrropuction. 


A method of computing approximately the coefficients of those terms of high 
order in the perturbative function which derive their importance from the pres- 
ence in them of a small divisor has been set forth by Porncar&. As applica- 
tions of the method to the problem of three bodies, Poincaré + has considered 
the case of coplanar orbits, one elliptic and one circular, and the case of 
coplanar elliptic orbits, both eccentricities being small. The latter case has 
been developed in detail by CocuLesco.{ Hamy$§ has discussed the cases of 
an elliptic enclosing a circular orbit and a circular enclosing an elliptic orbit, the 
mutual inclination being small, and the case of circular orbits with any inclina- 
tion. The latter of these was discussed anew by F&RAuD || together with the 
case of an elliptic enclosing a coplanar circular orbit, and the case of a circular 
enclosing a coplanar elliptic orbit. 

Owing to the high degree of the algebraic equations which define the singu- 
lar points upon which the approximate values of the coefficients depend, it is 
found necessary, in the general case, to resort to numerical examples. This 
paper is devoted to a discussion of these algebraic equations in certain cases 
more general than those so far considered, and to an extension of the criterion 
of admissibility as employed by Porncarté. The singular points in the most 
general case and in an important particular case are found and classified. The 
results in the cases studied by Porncar&é and CocuLESco are given first because 
of their usefulness in treating the more general problem. Before proceeding to 

* Presented to the Society (Princeton) September 13, 1909. 

{ PorncaRE. Les Méthodes Nouvelles de la Mécanique Céleste, vol. 1, Chap. V1. 

{CocuLEsco. Sur les expressions approchées des termes d’ordre élevé dans le développement de la 
fonction perturbatrice, Journal de Mathématiques, ser. 4, vol. 1 (1895). 

2Hamy. Sur les inegalités d’ordre élevé, Bulletin Astronomique, vol. 10 (1893). Sur le 
développement approché de la fonction perturbatrice, Journal de Mathématiques, ser. 4, vol. 
10 (1894) ; ser. 5, vol. 2 (1896). 

|| F&raup. Sur la valeur approché des coefficients d’ordre élevé dans le développement en series, 
Annales de l’Observatoire de Bordeaux, vol. 7. 
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this discussion, it is necessary to set forth briefly the principles of PorncaRé’s 
method, but in doing so the general equations are derived at once. 

It is required to develop a certain function F’? in terms of the mean anom- 
alies 7 and /’ in the form 


wo 
F’=> > hap i=V-—l, and m,, m, integers. 


m,=0 m,=0 


The coefficients A,,,,,,, are given by the integral 


1 —n—l 
1 


F(z, t) = 


where 
and 


The new variables, ¢ and z, have been introduced by the transformation 


1 
e = t’, = 


and the new constants by the relations 
m,=an+b, m,=en+d, 


where a, b, c, d are finite integers, small relative to n, and a and ¢ are rela- 
tively prime. Now if ®(z) is expanded as a Laurent series 


n=0 n=0 
and the result substituted in the expression defining A,,,,,, we get 


A mymge a, 


A theorem by Darsoux * can be employed to obtain approximate values of a, , 
and consequently of A.,,,,.,, for large values of n without computing the interven- 
ing coefficients. 

The initial contour of integration, |¢| = 1, may be deformed in any manner 
without affecting the value of the integral ®(z) provided no singular point of 
F(z, t), considered as a function of ¢ with z as a parameter, is encountered. 
By varying z continuously the singular points of /’(z, ¢) are displaced continu- 
ously. Singular points of ®(z) are those values of z for which two singular 


* DARBOUX. Sur V approximation des fonctions de trés-grands nombres, Journal de Mathé- 
matiques, ser. 3, vol. 4 (1878). 


d 
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points of #’(z, t), considered as a function of ¢, come into coincidence. If the 
points were initially, that is for |z| = 1, one inside and the other outside the 
contour |¢| = 1, then the singular point z of is admissible ; for, ®(z) 
ceasing to be holomorphic, the theorem of DarBovux is applicable. 
The principal part of the perturbative function * is 
F° 


where 
P’ = Bn)(£ — + 3(008 1) + % — (8,8 + 
L’, sin 1, uw, represent respectively the major axis, the eccentricity, the 
mean anomaly, the eccentric anomaly, and the longitude of perihelion of a planet 
B referred to the line of nodes and the center of gravity, D, of A and Bina 


system of three bodies A, B, C. L’*, sin ¢’, ete. are the elements of C 
referred to D and the lines of nodes. J is the inclination. &, 7 are defined by 


= cos u — sin + i cos dsin u, 

n = cos — sin + i cos sin wv, 
and &,, 7, are their conjugates. With this notation the codrdinates of B and C 
respectively are given by the real and imaginary parts of 

EL’ e, e’, 
respectively. 
The singular points of #’(z, ¢) with which we are concerned are the same as 


those of #’?. These are (a) the points for which w and wu’ as functions of / and 
l’ cease to be uniform, given by 


1 i 0 1 in ‘—0 
=1—sindecosu = = 1—sin cosu’ = 
Ou $ Ou 


and (b) the singularities of F’), given by 


19d 
P’ = 90. 
Making the substitutions 


wae", j=}(cosJ—1), 


=T, 1l+r=a, 


* This form of the perturbative function is due to FERAUD, loc. cit. 
Trans. Am. Math. Soc. 32 


72 
tan 5 =T, =o, 
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we get 
(1 
f= ox” ox 


Then the equations defining the singularities become 
(1) 
(2) (y—7)(1—yr) 
P = [o'y(x — Box(y—7')*] 
(3) x B,ow(1 — + joo'xy[y(a —7)(y — 
+ — 7)*(1— yr’ P— — — = 0. 


Singular points of the first species of ®(z) are given by the intersections of 
the curves (1), (2), (3). Singular points of the second species are those satisfy- 
ing the condition for double roots of (1), (2), or (3), any one of which may be 
written in the general form in the original variables 


(2s t) =Q. 
The condition for a multiple root of f(z, ¢) is 


of cou" of ao’ y? of 


Equations (1) and (2) will not give rise to points of the second species. 
Consider equation (3), which may be written 
P(x, BC+ joo'xyA, 


where 


(1 — yr’)? 
— B( —2r)(y—7'/, 
B=a'y(x—7) — Box(y—7'), 
C=o'y(1—27)'— Biox(1— 
The equation 0 P/0t = 0 becomes, after reducing by means of P = 0, 


R= —yT) 
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OA 
)| 4 | 


(4) 


OA 
+ jac'y (x A voy |- 0, 
where 


+ er)(y—1') + + )(w—7). 


§ 2. CopLaANaR OrsiTs. 


The equations P = 0 and 2 = 0, which determine all the singular points of 
the second species in the most general case, admit of important simplifications 
in certain particular cases. If we suppose the inclination zero, then 7 = 0, and 
we have 


(1) (3) B=0, (5) B=0, 
(2) (y—7r)(1—yr)=0, (4) C=0, (6) C’=0, 


where (3), (4) are the factors of the reduced equation P = 0, and (5), (6), 
those of 2 = 0. In this case the singularities are as follows: 


B=0, x=7,and y=7 intersectin UV, C=0,and y=7 intersect in V and W’, 


Cc=0, and intersectin U’, B=0,and intersect in V’ and W, 


1 
C=0,and x=T intersect in and Y’, x=T,and intersect in X, 
B=0,and intersect in Hand Y, x=-,and y=7 intersect in X’, 
T T 


B=0 and C = 0 intersect in B’, B, R, and PR, 
B=0 and = 0 intersect in C’, D, 7, and F, 
C = 0 and C’ = 0 intersect in 7”, D’, C, and F”. 
These singular points have been discussed by Porncaré and CocuLEsco. 


§ 3. THe GENERAL CASE. 


The main difficulty which arises in the general case is that of solving the 


algebraic equations. Even when 7, 7’, j are small, as is usually the case in 
nature, the approximations thus made possible do not so reduce the equations 
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that they admit of literal solutions. But the problem is not without interest 
from the theoretical standpoint, although the method it is necessary to employ is 
dependent upon numerical assumptions. In fact the results to be obtained in 
this manner have the same generality as those which would be furnished by a 
literal solution. 

Since for j = 0 the equations P = 0 and # = 0 reduce to 


B= 0, C=0, B=0,; C’ =0, 


it is easily seen that all the singular points for a small inclination will lie near 
those given in §2. The reducible curve P = 0 will lie near the curves B = 0 
and C = 0, but R = 0, being irreducible, will not lie near B’ = 0 or C’=0. 

Referring to equations (3), (4) of § 1, and recalling the fact that 8, 8, and 
also y, Y, are conjugate complex quantities, we can easily verify that 


A(x, wear 


Substituting in P and P above 1/x,1/y,1/i for x, y, i, multiplying through- 
out by «*y*, and changing signs in the latter function, we arrive at the identical 
expressions for P and #. Thus P = 0 and #& = 0 are self-reciprocal equa- 
tions in x, y, provided that 8, 8, and also y, y, interchange. This was to 
be expected, since the equations were originally in cos cos wu’, sinw, sin wv 
which involve x, y reciprocally. Consequently the discussion of the singular 
points will be very much abridged. 

The equation P = 0 is an algebraic equation of the sixth degree; the term 
with j as a factor is also of the sixth degree. The highest power of x, and also 
of y, is the fourth. The curve has thie axes as double asymptotes, and has also 
two oblique asymptotes. The equation 22 = 0 is of the seventh degree. As 
the highest power of x or of y is the fourth, the curve has three asymptotes 
parallel to each of the axes. Among these are the two axes. There is one 
oblique asymptote. 

We come to consider certain restrictions which may be put upon the elements. 
In the planetary theory we are concerned usually with the case of orbits having 
small eccentricities. We proceed to apply this restriction to equations P = 0 
and 2 = 0, and to note what simplifications take place. It was found in the 
simpler cases that the abscissas of half of the singularities were of the order 7 or 
7, and that the other half were of the order 1/7 or 1/7’. The same is true for 
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the ordinates of the singular points. If 7, 7’ are sufficiently small we have 
approximately 
o=1l+r=1, o=l+rTr=l1; 


and if in addition # and y are small, we have approximately 
(l—ary=1, 1—y7r =1, (l-yry=1. 
In the expression for A the coefficients of (y— 7’)? and (1 — y7 )* may be 
written respectively 


which become, omitting terms of the second order, — 8 and y¥, respectively. 
Making in P = 0 the approximations indicated and arranging in powers of y, 
we get 


P,= y[— Be(1+j)]+ + + 287'(1 +j)xr] 
The equation 2 = 0 likewise takes a reduced form. We have approximately 
OA 


oA 


Introducing these expressions in 2 = 0, simplifying as above, and arranging in 
powers of 7, we get 


R,=y[aB (a —7) 4+ che) + — 7) — 
+ y[—cB, x(a — 7°) — — — — — 2aB8,7'x( — 7)] 
+ cB, 7?) + 2aBB, (x + jey,7 2 
where 8’ = 8(1+,). 

Thus we have a system of equations, (1), (2), P, = 0, 2, = 0, for determining 
the singular points of ®(z) which lie near the origin. It is unnecessary to dis- 
cuss the geometry of the curves P, = 0, 2, = 9%. We know that in the region 
in which x, y are of the order +t, 7 these represent very closely the curves 
P=0,R=0. 


The equations for P,, 72, are of the forms 


Rh, = + BPy+hP=9, 
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where the coefficients a‘/), 8B are polynomials in x in which 7 is the highest and 
Eliminating y from these equations we obtain the 


i the lowest power involved. 


determinant 
a ) al? ) a?) a? ) 0 0 
(1) (2 (3) 2) 
0 a a, as 0 
1 2 
A. = 
1) ) 3) (3) 
By 0 l 0 0 
1) 1) (3 (3) 
0 v 0 0 
1) 1) (3) 3) 
0 0 0 0 0 B: 


This equation is of the fourteenth degree in «, the lowest term being of the 
second degree. Thus » = 0 is a double root; the other twelve roots give the 
singular points D, /, 7, C, B,, B,, R,, R,, U,, U,, U,, U,. 

The equation A, = 0 can be solved only when it is made numerical by giving 
particular values to the constants. The roots, which may be complex for general 
values of the coefficients, must first be separated and their values found by con- 
This procedure gives, indeed, the singular points 
A more fruitful method is 


tinued approximations. 
required but leads to no criterion of admissibility. 
to trace the original curves P, = 0, R, = 0 for certain values of the elements, 
and then, varying these elements, to note what changes take place in the singular 
points both as to their positions and as to their admissibility. 

The equation = 0 is homogeneous in c, a4; consequently a variation of 
the ratio c/a may be applied very easily. But in practice the initial choice of 
c/a will usually be sufficient. In fact, the important terms in the perturbative 
development are those corresponding to small values of 


where v, v' are the mean motions. We have approximately, for n very large. 


av+ev=0, 


from which 


The most important applications of the theory under discussion will be to the 
The particular values of the elements are 
Let L*, sin ¢, ete., be 


ease of Jupiter and a small planet. 
chosen in what follows with reference to this problem. 


the elements as defined above of the orbit of one of the small planets and let 
L”, sin ¢’, etc., be the elements of the orbit of Jupiter. 
8 = 2, then the successive reductions of —a/c become }, }; eo ete. This choice 


Assume the value 


c Vv 
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has the advantage of making all the singular points real in the case of zero incli- 
nation and coincident major axes, the necessary and sufficient condition being 
c/a <9. For the purpose of more accuracy in representing the curves P = 0, 
R =0 by means of P, = 0, &, = 0 we shall assume 7, 7’, 7 somewhat smaller 
than in the case in nature just suggested. We assume accordingly that 


a=—l, c= 2, |8| = 2, T= .03, 7=.01, j=—.0005. 
If the perihelia coincide at the node, we have 
B=8,=7=% = 2: 
if at 7/2 from the node, 
On introducing these values, the equations ?, = 0, 2, = 0 become 
P, [—2a] + .0009] + y[ —2a*+ .04a° — .0020x = .00102] 
+ .0004a? = 0. 

R, = y'[ 2x + .06 ]+ — .0086 | + y[ — — + .0068x + .0030x 

— .000024 = .000030] + .04x* — .00082° — .0000162 = .0000202 = 0, 


the double sign arising through y = y, = + 2. 

In the table on the following page are the equations resulting from particular 
values of x near the origin. The roots of the equation are multiplied by 100 and 
the corresponding coefficients and roots are set down. In the first column the 
values of » are given; the two single and two double columns following give the 
coefficients of »*, y*, y and the constant term of the reduced equation. The three 
triple columns give the roots of the cubic in y. They represent respectively the 
roots for 7=0; j = — .0005, y= y, = 2: j= — 0005, y= = —2. 


§4. THe CRITERION OF ADMISSIBILITY. 


It has been stated in $1 that for a singular point of ®(z) to be admissible 
two singular points of /’(z, ¢) must have come into coincidence from opposite 
sides of the contour |¢|=1. But /’(z,¢) is a multiform function of z, ¢ of 
an infinite number of determinations. For /’(z, ¢) contains the expression P 
uniform in a, ¥, which are in turn multiform in z, ¢, as is seen from the relations 


The discussion for admissibility must then be made for each determination of 
F(z, t), that is, upon each sheet of the corresponding Riemann surface. 
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For the particular case t’ = 0 Porncar& has overcome the difficulty. Using 
this value of 7’ in z above and solving for y, we get 


1 


in which it is seen that y is uniform in 2’”, #'”, i. e., in 2’, x’ if we put 


Also P(x, y) is uniform in z’, x’, and in these variables /’(z, ¢) has only two 
determinations, namely those of the radical P-!. This is effectively a reduction 
to one sheet, as these two determinations are equal in numerical value and oppo- 
site in sign. In the new variables the initial contour |¢| = 1 becomes |#’ | = 1. 
This method has been extended by CocuLesco to the case of small values of 7’, 
where as a first approximation the eccentric and mean anomalies are considered 
equal, 

In general it can be shown by expanding 2’ as a power series in 7’, and neg- 
lecting terms above the nth power (as is permissible except for values of y 
of the same order as 7’, such as occur along certain branches of the curve 
P=0) that F(z, ¢) becomes a function of multiplicity 2” in the variables 


z',”'. The geometrical meaning of this transformation as applied to the con- 


tour |¢| =1 is interesting, but the investigation is unnecessary to the present 
discussion. For it will be shown that the admissibility of points for the general 
case of tT’ + 0 can be determined from the discussion of the case of tT’ = 0. 

In fact, Porycaré& has shown that, in general, a singular point of ®(z) can- 
not change its character as regards admissibility for any variation of the ele- 
ments. The exceptional cases, as we shall see, are those for which there arises a 
change in the combinations of the singularities of /’(z, ¢) which form the singu- 
larities of ®(z). We proceed to establish a theorem simpler in form than the 
original theorem by Porncar&, but embodying its essential point. This theorem 
may be enunciated thus : 


I. For the values |z|=1 the singular points of F(z, t) cannot cross the 
circle |t} = 1. 

Having seen how the circle |¢|= 1 can be replaced by either |x| =1 or 
|a’ | = 1, we are at liberty to make the discussion in the plane of x. Consider 


the relation 


If we let z be a parameter, this relation is an implicit multiform function of the 
two complex variables x and y. We have seen that it has an infinite number 


of determinations. There corresponds then to the function z(x”, y) a Riemann 


Yor 9 
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surface S with an infinite number of sheets. The surface varies when the par- 


ameter z is varied, and for |z| = 1 becomes the surface S,. Upon this surface 


trace the circle C,, defined by 


| 
x | 


| y|=1. 

This is possible, for if « be given a value of modulus unity then we can choose 
a value of y having the modulus unity in such a manner that z may have any 


value whatever of modulus unity. On the other hand if |x| =|y|=1 and 


x = y, then there corresponds a unique point on the circle of unit radius in the 
z-plane. 

Now to each point on the Riemann surface S there correspond definite values 
of x and y, and consequently definite values of the eccentric anomalies w and w’. 
The points on the circles 


|y| = 1 

give real values of the eccentric anomalies, i. e., a real point on each of the 
orbits. But in this case #’(z,¢t) has a singularity only when the distance 
between the two bodies is zero, that is, when 


P=0; 
or, when the distance of one of the planets from the sun is zero, that is, when 


1—sindeosu=0, 
or 
1 — sin ¢ cosa’ = 0. 


Excluding then the case of intersecting orbits and also the limiting case of the 
problem of two bodies where the orbit reduces to a straight line, it is impossible 
that for any variation of the elements a point on the circle C’,, shall become a 
singular point of /’(z,¢). Therefore a singular point of #’(z, ¢) can not cross 
the circle |¢| = 1 for values |z| = 1. 

A singular point of ®(z) has been defined as being admissible when it is 
formed by the coming into coincidence of two singularities of /’(z, ¢) from the 
interior and exterior respectively of the initial contour |t]=1. The theorem 
just proved may be applied to the question of the permanence in the character 
of admissibility of the singular points of ®(z). To distinguish among these 
singular points certain definitions are useful. We assume that the elements of 
the orbits, with the exception of the anomalies, remain constant. Let A, B 
be two singular points of ®(z) at which 2’ has the values z|, z, and assume 


|z|<|2;|<1. When 2’ is made to vary starting from the point A, then the 


corresponding points in the plane of x’ separate and move subject to the condi- 


tions defining singular points of #(z,¢). If in the plane of a’ there is a 
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means of direct passage from A to B along a line of continuously increasing 
modulus of z, then B is subordinate to A.* The means of passage is the locus 
furnished by one of the equations (1), (2), P = 0. 

If the point z’ encounters no other singular points in its course, and if the 
points A and B are only double points, then A gives three terminal values in 
the plane of 2’, one proceeding from A directly and two proceeding from the 
branch separating at B. ha 

We define a singular point as primary if it has no admissible subdrdinate 
point. It is secondary if it has an admissible subordinate voint. A primary 
point may be either admissible or inadmissible; a secondary point is pseudo- 
admissible, for, although there correspond to it terminal points both within and 


without the circle || = 1, there can be no change of elements which will bring 
it upon the circle of convergence |z’| = 7’ without also causing it to pass through 


the circle of absolute value of its subordinate point, and the subordinate point 
itself determines the limit of convergence of the Laurent series. 

With these definitions we proceed to the discussion of admissibility, taking 
up the different cases that may arise. It is requisite first to establish the 
theorem : 

Il. A primary singular point of ®(z) can not change its character as 
regards admissibility for any variation of the elements, and this character can 
be determined by varying 2 along any path whatever of increasing modulus. 

Consider the method of formation of the singular points of ®(z). The vari- 
ables a and y satisfy certain equations of the series (1), (2), P= 0, 2 = 0. 
The intersections furnished by these equations give singular points of ®(z) of 
the first and second species, and to such a point there corresponds a definite 
value of z’, that is, a point in the plane of z’. Instead of finding the intersec- 
tions indicated and then computing the values of z’, it is simpler to introduce 
certain auxiliary curves to take the place of (1), (2), P = 0, R = 0, whose 
intersections give the required singular points of ®(z). Let 


a(x,y)=90 


be any one of the curves (1), (2), P=0, R=0. We have 
2(2,y)=r. 


When »’ varies in a given manner, the points of intersection of a= 0 and z =r 
describe certain curves; in particular, if 2’ describes the circle 


, 
| 


=T, 


|z 


the singular points arising from the intersections of the curves will cause 2 to 


* POINCARE, loc. cit., p. 308. 
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describe certain 8-curves in its plane. The curves a = 0 may be replaced by 
the corresponding 8-curves. Mere intersections of the 8-curves are not singular 
points of ®(z); the tracing points must coincide for some value of 2’. 

For the initial circle |z’| = 1 the §-curves in the plane of x’ do not intersect 


the circle |x’| = 1, by theorem I. Now let the circle |z’| = r’ contract, keeping 


the origin fixed. The 8-curves sweep over certain portions of their plane. The 
significance of this is that corresponding to any continuous path whatever from 
= yr 


, 


any point on the outer circle |z’| = 1 to any point on the inner circle 


there are certain paths lying entirely within the region swept over by the A- 
curves. Suppose that when the circle |z'| = 7’ is reached two tracing points of 
the 8-curves coincide. Corresponding to this point in the a’-plane is a point on 
the circle |z'| = 7’, namely one of the singular points of ®(z). Now reverse 
the variation of z’, and starting with the singular point let z’ go to a value 
|2'| = 1 along any path whatever consistent with the continuous increase of the 
radius of the circle on which it is found. The points in the a’-plane accordingly 
trace definite loci in the swept-over area and rest finally on the initial S-curves 


corresponding to |z’'|= 1. These curves being entirely within or entirely with- 
out the circle |x| = 1, the terminal values of 2’ are definitely within or without 
the circle || = 1 for any z-path to the circle |z’| = 1; whence we conclude 


that the point can not change its character as regards admissibility for any 
change of the elements, for this would necessitate that some of the points of the 
8-curve should cross the contour 


which we have shown to be impossible. 

If in a series of singular points, connected by relations of subordination as 
defined previously, the one first encountered is inadmissible, the cirele |z' | = r 
may be contracted over this singular point, since a deformation of the contour 
in the plane of a’ is not rendered impossible. This may be repeated until an 
admissible point is reached. It is evident that the theorem applies to this point 
as well as to all the preceding points of the series, that is, to all the primary 
points of the problem. 

Consider now all the points both primary and secondary in the z-plane con- 
nected in any manner whatever as required by the definition of subordination 
and as permitted by the method of varying z. If the elements be varied con- 
tinuously these singular points likewise vary continuously. Their curves of con- 
nection must be deformed if necessary to correspond to the displacements. It 
is proposed to account for the behavior and the consequent character of the 
points for all possible changes on their relative positions. 

Let A represent a typical singular point with B as an immediate subordinate 


= |y| = 1, 
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point, and consider any particular series of connected points. Let B pass 
through A. Then: 

(1) If both were originally inadmissible they remain so after the change, for 
their combined terminations in the plane are still on the same side of the unit 
circle. 

(2) If A were originally admissible and PB inadmissible, the passage of B 
through A leaves A admissible and renders PB admissible and secondary, i. e., 
pseudo-admissible. 

(3) There is no case of B admissible and A inadmissible. 

(4) If B is admissible, then A is pseudo-admissible, and the passage of B 
through A gives the following result : 

(a) If A were intrinsically admissible, that is, if its admissibility did not 
depend upon a branch of B not common to A and B, then A remains admis- 
sible and 2 becomes pseudo-admissible. 

(6) If A were not intrinsically admissible, then A ceases to be admissible, 
and B remains admissible. 

This discussion covers all cases of interaction of primary singular points. 
This interaction produces certain secondary points. If we consider the passage 
of a secondary admissibility through an immediate subordinate primary admis- 
sibility, we shall have one of the above four cases, consequently the converse 
results hold. In fact, these converse relations are contained in the cases as 
already stated. We have: (1) is its own converse, (4) (a) is its own converse, 
(3) does not exist, (2) is the converse of (4) (0), 

It has been assumed that B actually passed through A, but the results would 
be the same if it passed by it simply, i. e., passed through it in absolute value. 

It should be remarked here that the variation of z from a secondary point to 
to = 1 is not entirely arbitrary. For if all subordinate points are not 
taken account of, it may easily happen that we get an apparent inadmissibility 
when the point is actually admissible as found by the path through the sub- 


vA 


ordinate point. 

This completes the theory of admissibility. Before proceeding to the discussion 
of the points found in § 3 it should be remarked that there are for each point 
in addition ec — 1 others of the same modulus and differing in amplitude by 
the ares ka/c. Portncar& showed that these points are all admissible or inad- 
missible according as one of their number is admissible or inadmissible, and 
this holds true in the general case. 


§ 5. Discussion OF THE SINGULAR POINTS. 


Having seen that the singular points which arise in ease of particular elements 
can not, as a rule, change their character as regards admissibility for a change 


of elements, it is advantageous to make use of the simpler cases to determine 
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what points are admissible and to follow these points when certain of the ele- 
ments are made to vary continuously, in order to determine in the general case 
the points that are admissible. 

Accordingly we make direct use of the results of Porncaré and CocuLEsco 
already cited. Let us designate three cases as follows: Case I, 7 = 0; Case 
II, 7 + 0, perihelia coincide at a node; Case III, 7 + 0, perihelia coincide at 
ninety degrees from a node. In the figures 1, 2, 3 the light and heavy lines 
are the curves in the case 7 = 0, 7 + 0 respectively, also the full lines are the 
curves P,=0,2=7,y=7, and the broken lines are curves 2, =0. The 


Fig. 1. 


positions of the small ares show where, with respect to the circle |a’| = 1, the 


z|=1. 


points on the curves ?, =0,2=7, y=T lie for which 


Case 1. j=0, ©=@. This is the case discussed by Porncarfé and 
CocuLesco. The singular points have already been classified in § 2. Of the 
points of the first species B, 2, #, W, X, ¥ are inadmissible and U, V are 
admissible. Of the points of the second species, C’, /’ are inadmissible and D, 
T are admissible. 

The variation of z is as follows: On the branches 8,2,, z decreases to C, 


which in this case coincides with 2, is a minimum at C, increases from C’ to O, 
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changes at O from |z| = 0 to |z 
the branches 8,8,, %,%,,2 decreases to D, is a minimum at J), increases through 
B to P, |z| =o at P, = 0 at Q, increases to F’, which in this case coin- 
cides with /2, is a maximum at /’, decreases from F’ through V. On the 
branch 8,4,, 2 increases to 7’, is a maximum at 7’, decreases from 7’ to O, 
z| = 0 to |z| = 00, decreases from O to U, and increases 


2|= 


= 0, and increases from O through #’. On 


» 


changes at O from 
from U. Along x =7T the value of z increases with y positive, and along y= 7’ 
it increases with a positive. 

Case I]. j4+0,@=0'=0. The curves in this case are given in Figures 


1, 2, which are plotted with the values given in Table I. The branches a,, a,. 


2. 


ete., are the curves of the Case I, just considered ; the branches a, a), ete., are 
for present consideration. The curves ?,=0,2=7, y=T are represented 
in full lines, and the curves 22, = 0 by dotted lines. 

Consider the changes which the original curves undergo when the inclination 
is made to vary from zero, in other words when j varies from zero to a small 
negative quantity. The segments a, #2 and 2O vary continuously to form the 
branch The segments a, 2, RU, Ua, form the branch «,U,a;. The 


segments OU, Ua, form the branch OU,a,. The straight lines y= 


‘ 
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are of course unchanged. The branches a,a, and a,a, become the branches 
a, a, and a, a respectively. 

The singular points for j sufficiently small will lie very near those of Case I. 
The point £ is in reality a combination of what would be in general the points 
Fand &. This coincidence is due to the choice of constants in the particular 
numerical case computed. A change of the ratio c/a is sufficient to separate 
the points /’ and # without changing the curve P = 0 on which they lie. But 
these points can be separated otherwise. The branch a, 2U is tangent to the 
branch a, RU; this gives two coincident points of intersection at 2. The 
branch a, 2 U also cuts the branch a, RO at F in a point of the second species 
R,. Thus F is composed of the points /’, 2,, 2,, all of the second species. 
When j varies in the manner indicated the point /’ remains rea] and takes a 
position /’,; the two points 2,, 2, become conjugate imaginaries. 

The point U in the simple case is composed of nine points. These arise as 
follows: each of the two branches of 2 = 0 which pass through U cuts each of 
the branches of P = 0 at the same point, giving rise to four points U,, U,, U;, U, 
of the second species; each of the straight lines x = 7, y= 7 cuts each of the 
branches of ? = 0 in the same point, giving rise to four coincident points of the 
first species, U,, U,, U,, U,; the point U is also the intersection of x = 7, 
y=t. Now vary j as before; the nine points all separate. U, and U, of the 
second species remain real. UU, remains real. U, is imaginary in the figure, 
but it is evident that for sufficiently small values of j it is real. In fact it 
leaves U as a real point and different from the point /, these two points 
approaching each other along the line x =7. For a certain value of j the 
points # and U, coincide, and when / is given still larger values they become 
conjugate imaginaries. The points U,, U, are conjugate imaginaries; the 
points U_, U, are also conjugate imaginaries. 

The remaining points V, of the figure, and W, Y, X, not represented in the 
figure, take positions near their respective original positions. 

Discussion. — By virtue of the theorem established in § 4 it is not necessary 
to discuss the points R,, V7, W, X, Y, since we have already 
determined the nature of the corresponding points in the simpler case. We 
have only to consider the nine points corresponding to U. 

The equation P = 0 is of the form 


P= ST + jA. 


We proceed to find the form of this curve for the values of x lying between 
some value x of the order 7, and the valuea#=1. The product S7 considered 
as a polynomial in y is of the fourth degree, while A is a polynomial of the 
third degree in y. By taking j sufficiently small we can make the curve 


ST +jA=0 
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lie as close as we choose to the curve 
ST=0 


throughout the region under consideration; for the coefficient of the highest 
power in y remains unchanged, and all the coefficients of A are either finite or 
small. 

Then for a given value of x the corresponding values of z' along similar 
branches in the two cases will be approximately equal. In the case of j7 = 0 
the terminal values of x’, when z’ described the straight line of constant ampli- 
tude, were definitely within or without the circle |z’|=1. It is possible in the 
present case to choose j small enough to make the terminal value of x’, relative 
to the circle |x| = 1, correspond exactly to those in Case I. The character of 
the real points U,, U,, U,, U, is now evident. 

U, is formed by the intersection of x = 7 with a,U,a,. The two points on 
separation follow their respective branches in the direction of y increasing. They 
reunite in #. But £ is inadmissible, therefore U, is inadmissible. 

U,, is inadmissible. The two points which separate are found to follow 
U,O and x = T in the direction of y decreasing. The terminal values of x are 
both without the circle |z’| = 1. 

U, gives 2’ a minimum. When 2 varies along the positive part of the real 
axis of its plane toward the value z’ = 1, the two points x’ on separation move 
along U,O and Ua) respectively. The latter encounters in its course the 
point U,. All the terminations for UV; have been found to lie within |x| =1. 
The same is true for the terminations on the branch U,a,, therefore U, is 
inadmissible. 

U, is admissible. Since at this point z’ is a minimum the values of 2’ on 
separation follow the branches U,a; and U,a;. When 2’ attains the value 
unity we have for the first point « < 1 and for the second x > 1. 

The points U,, U,, U,, U, are left for consideration. These points are all 
imaginary but their amplitudes are very small when j is put small. Correspond- 
ing to each of these points the value of z’ is imaginary with a small amplitude. 
Instead of varying z’ along a line of constant amplitude, as was done in the case 
of the real points, let us vary z’ along the straight line limited by the points z; 
and z’ = 1, where z) is the value of 2’ at one of the singular points. The loci 
of the singular points which separate when z’ is made to vary from the value z; 
are certain curves satisfying the relations 


P(x,y)=9, 


for points of the second species U,, U,; or 


z(x,y)=6, 


Trans. Am. Math. Soc. 33 
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for points of the second species, U,, U,, where ¢ is a parameter varying along 
the straight line limited by 2; and z =1. In the case of j7 = 0 this variation 
of € is along the real axis, and consequently the singular points on separation 
follow P(x, y) =0 in the real plane. For j very small the singular points will 
describe curves lying not in the real plane, but very near the real part of P = 0. 
The approximation becomes closer as € approaches = 1, for which value the 
singular points have positions on the real part of P = 0 not far from those 
terminal positions found for j = 0. 

Consider the singular points of the second species, U,, U,, formed by the 
intersections of the curves P =0 and #=0 in the neighborhood of the 
point U7. Such an intersection implies that z’ is a maximum or a minimum. 
In each case we know it to be a minimum. Then the four points a’ after 
separating, two from U, and two from U,, follow each of them one of the four 
segments of P = 0, which approach coincidence with UO, UR, Ua,, Uz, 
respectively. 

An examination of the order of interchange of the branches near U is neces- 
sary. But this discussion is practically obviated through @ priori considera- 
tions. In fact, U7 represents an ordinary node of the function ? = 0 in the 
case of j = 0. When j varies the interchange may take place in one of four 
different ways. The double point may continue to exist near U, or the two 
branches may cease to intersect and the continuous branches may lie near the 
original curve in any of three ways as follows: a;’UR and a,’UO, a;'Ua"’ and 
RUO, a and It can easily be shown that if the node continued 
to exist for variations of j, there would be one or three admissible singular 
points according as one of the two points of the second species is or is not 
on the segment U\,<;’. But it can be shown that for general values of 7 
there is no node U,,. Consider the three other possibilities. By virtue of the 
continuity of the function P we see that each of the four segments which 
the singular points describe after separation is involved once and only once. 
We know that of the four terminal values of « three are within the circle 
|| = 1 and one (on the branch 2;’) is-without this circle. Thus of the two 
conjugate imaginary points, U/,, U,, one is admissible and one is inadmissible. 
This distinction will be considered later. 

The points U,, U, which are the imaginary intersections of y=7' and P =0 
in the neighborhood of U7 become real when @ = @ = 7/2, and are reserved 


for discussion under Case III. 

The point which is the intersection of « = 7 and y= 7’, is or is not 
admissible according as one of the points U_,, U,, supposed subordinate, is or is 
not admissible. This illustrates a distinction made in the foregoing section. 
In fact, if the points on separation were permitted to follow the line y = 7’ only 
in the real plane, we should be led to conclude that U is necessarily inadmissible. 
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The curves as they exist in the third quadrant are represented by Fig. 2. 
The branches £,, 8,, ete., are the curves of Case I, while 8/, 8), etc., are the 
curves in the present case. From the figure it is seen that as j varies from zero, 
the line 8, TOBS, becomes 8;, and 8, DBB, becomes 8,8). The branch 
8, DTOBS, becomes 8,08'. At B three points coincide. These separate 
when j is made to vary, the point C remaining real, B, and B, becoming con- 
jugate imaginaries. The points D, 7’, remain real and their abscissas increase 
in absolute value. 


Fig. 3. 


The points B,, B,, C are inadmissible and the points D, 7 are admissible, 
as was found in Case I. 

Case HI. j + 0,@=@ =7/2. In this case the curves as they exist in 
the first quadrant are represented by Fig. 3. They are plotted with the partic- 
ular values given in Table I. The branches a,, «,, etc., are the curves in the 
case 7 = 0; the branches |, a), ete., are for present consideration. 

When j varies from zero to the small value used in the table, the original 
curves vary as follows: a, Ra, becomes a,a,; O/UO remains a closed curve 
and recedes toward the point O; a,Ua, becomes 2/a/; the branches a, 2 Ua, 
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and a,Ua, become respectively the branches a} Ua, and a,a,; finally, the 
branch a,O passes through a loop form, then to the position Oa. The order 
of the separation of the branches here is seen to be opposite to that of Case II. 

The singular point at 2 breaks up into the points /’, R,, R,. The points 
F’, R, are conjugate imaginaries. /?,, the point lying below /?, is real at first, 
but for the value of j employed in the figure it has become imaginary. When 
the loop formed by the segments OR UO begins to recede toward the origin, 
the two points of intersection of this loop with a line RU, little different from 
the lower of the dotted lines 2U, will be real and different. Designate by U, 
that one of the points near 7. The value of j continuing its variation, the 
points U, and 22, come into coincidence and then separate to become conjugate 
imaginaries. Consider the variation of z’ on this loop when U, and £, are 
real. The loop is cut by the line y= 7’ at the points V and U,. The varia- 
tion of z’ is as follows: zero at O; OVAR,, increases; maximum, at 2,; 
R,U,, decreases; minimum at U,; U,U,O, increases; infinity at O. Now 
let the loop decrease, the maximum point #2, comes to coincide with the mini- 
mum point U,, then z’ continuously increases around the loop OV,U,O. The 
loop continuing to decrease, V, and U, approach each other, coincide, and then 
become conjugate imaginaries. The curve a; Ua; and the line y= 7’ cut the 
curve P = O at approximately the same point in the figure. The straight line 
x2 = T cuts the curve P = O in two imaginary points U,, U,. 

Discussion. — By using the results of the two preceding cases the character 
of the singular points in the present case is easily determined. 

F is inadmissible. 

E is inadmissible. 

f, has F' and £ as subordinate points and is itself inadmissible. 

f, is a maximum point of z’, and when 2’ increases the two points become 
conjugate imaginaries at first and then coincide at #,. Thus R, is inadmissible. 

U,, the intersection of a{ a, and Ua,, is admissible. For, after separation, 
the two values of 2’ follow the branches a; and a) respectively, the terminal on 
the former branch being without and that on the latter branch being within the 
circle |x’ | =1. 

U,, subordinate to U, for sufficiently small values of 7, is admissible. The 
tangent to the curve Ua; at U is negative thus giving U, a position above U,. 
The terminal value of x’ in the directions of x increasing and a, are respectively 
within and without the circles |x’ | =1. 

U has U, as a subordinate point and is consequently admissible. 

U, is admissible ; of its two branches one follows the curve U,, O and the 
other passes through U and U, which are admissible. 

U, is admissible; it has #, and U, as subordinate points, of which the latter 
is admissible. 
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y 


V, is admissible ; it has at first 2, and U_ as subordinate points, and later U : 


alone. 

U, and U,, the intersections of « = 7 with P = O are admissible and inad- 
missible respectively. For, of the four branches UR, UO, Ua,, Uz,, only the 
latter gives a terminal value of x’ greater than unity. With the argument used 
for U, and U, of Case II, the conclusion is evident. The discussion for the 
third quadrant is not difficult and may be omitted. 


$6. Tue Case or 7 =0,7+0. ConcLusion. 


Before passing to conclusions relative to the general case, it is interesting to 
consider what simplifications arise on taking the eccentricity of the orbit of the 
second planet zero, the eccentricity of the orbit of the first being small, and the 
mutual inclination small. 

Neglecting in the equation P = 0 the terms in 7 and arranging in powers of 
y, we have 


+j)]= 0. 


An approximate expression for this curve in the neighborhood of the origin can 
be obtained in a very much reduced form by neglecting all but terms of rela- 
tively low order in the various coefficients. Thus considering « of the same 
order as 7, and neglecting all but terms of the second order in the coefficients, 
there results the equation 


P, = + + + = 0. 
Likewise the equation 2 = 0 becomes 
R, = —7)(1 + Box(1 +j)] 
+ [—a8,(x—1)— cB, — 7?) — cy,ox + ay,(~—T)}] =9. 


The eliminant A, = 0 of P, = 0 and #, = 0 may be found as in $3. It has for 
its roots the abscissas of points corresponding to U, F’, R,, R,, D, C, B,, B,. 
Equations P, = 0, #, = 0 may be written in the expanded forms 


2] + + + + [850° + + = 0, 
[tye +t, J+y[t.2+ toy] + [ + + + tro = 0, 


t,,, t;, are so designated because they contain all the terms in j, and 


where s,,, 
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to these terms the double sign is attached. Then A, becomes 
[ + + + utat+-.- ] 
+ [t,,2° + + + 
where the w, v are functions of the s, ¢ as indicated. Or again, 


Assuming the approximate values of the roots as found by plotting, they may be 
corrected by applying the well known formula 


A,(,) 


which is true for all simple roots, real or imaginary. The work may be reduced 
in the two cases under consideration by finding the four real roots first and then 
reducing the order of the equations. 

In the two cases of o = @ =0 and ®=@' = 7/2, the roots of A, may be 
calculated very rapidly by the above formula when j is varied little by little, 
t,,, and w,, v, in an additive manner simply. 


since changes in j will affect s,,, 

Since the choice of 8, 7 can be made once for all, and since j can be made to 
assume any assigned value in the manner indicated, there remains for considera- 
tion the variation of the longitudes of the perihelia. Having taken 7’ =0, we 


have ® = @ and 
72 


B=B,= 
When G is taken different from zero, imaginary terms affect the expressions A,. 
The roots of A, = 0 become imaginary. They may be computed by the above 


formula as @ varies from zero. The computation of «x in the functions 
A,(x,), A°.x,) is best effected by the use of trigonometric tables by writing 


= |x, |"e™? = |x, |"(cos nO + isin 


From the foregoing case the means of extension to the most general case is 
now evident. We have seen that it is necessary to take account of the analytic 
possibility that a singular point may become admissible or cease to be admissible, 
but slight attention to the geometry of the variation of the singular points as 
the elements vary will suffice to guard against error. Probably the selection of 
the primary admissible points, as well as their tracing in the x-plane, with the 
application of the principle of the limits of the roots of an equation, will answer 
in most of the practical applications. Approximations may usually be made in 
the manner employed here, and where necessary recourse may be had to the 
original equations of a higher degree. 
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Only points within the circle |z| = 1 have been considered in the foregoing 
discussion. After we have determined the coordinates ., y of the singular point 
which gives z nearest the circle, the point outside the circle || = 1 and nearest 
it is found as follows. We have 

A(x, i)= XA ( —i)=0, 
where Y isa power of .. and A is the equation in « determining singular points. 
This shows that 1/2 is also a root of Aif —iis put fori. But this is effec- 
tively a change of elements, namely an interchange of 8 with 8, and of y with y¥,. 


Make this change, then 


1 
— xa( i)=0. 


Thus we see that by putting — i for i the knowledge of the inner point which 
conforms to the demands of the problem enables us to find the required outer 


point. 
The simplest relations from which to begin the variation of the perihelia are 


with & coincident or opposite, which gives 
L” 
8(=8,)=+7(=%) == L?° 


The variation of c and a has not been considered. These constants enter the 
equation 22 = 0 only. This equation is homogeneous in c, «, so that if it is 
decided to vary the ratio c/a after a singular point is found, a factor may be 
entered readily in each term. Otherwise the problem may be begun anew for 
each choice of c/a, c and a being small integers. 
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NOTES AND ERRATA 
VOLUME 8 


G. A. MILLER: Generalization of the groups of genus Zero. 

The theorems in italics on pages 12 and 13 should read as follows: There 
ave exactly six non-abelian groups whose two generators s,, 8, satisfy the equa. 
tions = ( 8, = 1 They area group of order 60.32 = 1920 and the 
direct products of the icosahedron group and the cyclic group of order 2°, 
4. 

There ave exactly sie non-abelian groups whose two generators 8,, 8, satisfy 
the equations = 83, (s,s, =1. They are the icosahedron group, G,,,, and 
the direct products of these groups and the cyclic groups of orders 5 and 25. 

In the published theorem corresponding to the former of these two, the group 
of order 1920 was omitted owing to the incorrect conclusion that the order of ¢, 
could not be 160. The error which is corrected in the second of these theorems 
resulted from a mistake in the formula below the middle of page 12, which 
should read 


In the published formula the last exponent was erroneously given as —12”. It 
should also be added that the theorems of this paper relate to non-abelian groups 
only, as the cases where the operators are commutative are so evident as to appear 
almost trivial. The restriction to non-abelian groups, however, should have 


heen stated distinctly, but this was not done. 
VOLUME 10 
Ie. KASNER: Natural s of trapectorie 


P. 262, i. 17. For FouRIER’s read FERMAT'S. 
244.46 Np. the ‘a. 


) 


